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topologies and its invariances allow to classify all CA rules into equivalence classes. Complexity 
in ID systems is then shown to emerge from the weak symmetry breaking of the addition 
' I— ' modulo an integer number p. The latter symmetry is possessed by certain rules that produce 

^^ Pascal simplices in their time evolution. These results elucidate Wolfram's classification of CA 
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Summary of the main results and conclusions 

Cellular automata (CA) constitute paradigmatic models of complexity in nature, from 
snowflakes, patterns in mollusc seashells and spiral waves in the Belousov-Zhabotinsky reaction 
to neural networks and the fundamental physical reality [TH3]. CA serve as models for complex 
natural systems made of large numbers of identical parts fW6|. They consist of a discrete lattice 
of sites, with a finite set of possible values each, although they even allow the dynamics of con- 
tinuum systems to be accurately described [7]. The site values evolve synchronously in discrete 
time steps according to identical rules, being determined by the previous values on the sites of 
their neighborhood. The concept of CA dates back over half a century to the efforts of John 
von Neumann to design self-replicating artificial systems capable of universal computation [8J 
(see also 0). It was however in the early eighties, when Wolfram published a series of ground- 
breaking papers on the subject [3H6l [TOl [TT] . that CA received wide attention from the scientific 
community. Wolfram classified C A into four classes of increasing complexity [4] . For a random 
initial condition, he found through computer experiments that a CA can evolve into a single 
homogeneous state (Class 1), a set of separated simple stable or periodic structures (Class 2), 
a chaotic, aperiodic or nested pattern (Class 3) and complex, localized structures, some times 
long-lived (Class 4), see FigfTl Despite some monographs and further important work [12H22J 
and tantalizing hints |14j . no theory provided neither a satisfactory explanation for Wolfram's 
observations nor analytical expressions valid for arbitrary neighborhood sizes, number of site 
values and dimensions. Mathematical tools for these discrete systems with the practical value 
that partial differential equations have for continuum systems [IT] were absent. This has been 
a major drawback precluding significant progress: although quite a tour de force, no general 



conclusions could be sharply drawn other than the computational ones [H [T^ [T5] . It has indeed 
long been thought that no simple mathematical expressions could be given for most CA and a 
strong emphasis has been put on computational aspects instead of standard mathematics [1]. 
This preprint tackles this problem providing the missing mathematical tool for CA. A universal 
map encompassing all ID deterministic first-order in time CA is presented here. This is to be 
considered as the discrete counterpart of partial differential equations in continuum systems. 
The map can be easily extended to an arbitrary number of dimensions and topologies. The 
advantages of having such a mathematical expression become soon apparent: symmetry argu- 
ments applied to the map allow to classify all dynamical CA rules into equivalence classes and 
a most surprising result is that a theorem can be proved which establishes how a CA rule is 
constructed in terms of rules of lower range, relating its behavior to them. This theorem and 
its consequences allow the rationale behind Wolfram's classification sketched above to be estab- 
lished. The crucial result is that the most complex rules (those with Class 4 behavior) can be 
found with a simple prescription, starting from rules possessing the symmetry upon addition 
modulo an integer number p, and weakly breaking this symmetry through the additional degree 
of freedom involved in constructing a rule of higher range. I illustrate how this mechanism is 
the origin of complexity in ID systems. 

I consider a general rule whose action is described by giving a table of configurations in- 
dexed by an integer n so that to each configuratioun corresponds an outputn- The rule is 
deterministic and each of these configurations is exclusive, i.e., a given input cannot be simulta- 
neously equal to two different configuratioUn within a certain tolerance = {configuratiounj^i — 
configurationn)/2. The output of the rule is then given by 

output = (1) 

y outputn ■ B{configurationn — input, tolerance) 

ndtable 

where the function B{x, e) returns one if — |e| < x < \e\ and zero otherwise and, thus, coincides 

with the boxcar function 

^, ■. ^ f X + e X — e\ 

B(X, e) = - ; r - 1 r (2) 

^ ' ^ 2\x + e x-e J ^ ^ 
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FIG. 1: Spatiotemporal evolution of one-dimensional CA rules representative of each of the 
Wolfram classes calculated from Eq. g. Rule ^254^ (Class 1), rule ^232^ (Class 2), rule ^30^ 
(Class 3) and rule ^1102 (Class 4). Time flows from top to bottom. 



(see Fig. SI in the Supporting Information |23]). 

I show now how CA rules fit always in this structure. I focus first on a ID ring containing 
a total number of Ng sites. An input is given as initial condition in the form of a vector 
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\ . Each of the Xq is an integer in the range through p — 1 where superindex i 



specifies the position of the site on the ID ring. At each t the vector X( = {x\, ...,Xf. ") specifies 
the state of the CA. Inputs and outputs from the rule are integers on the interval [0,p — 1]. 
Periodic boundary conditions are considered so that x^ " = x} and x^ = x^". Let xl_^i be 
taken to denote the value of site i at time step t + 1 . Formally, its dependence on the values at 
the previous time step is given through the mapping xl_^_i = (t){xl~^, x\~'^~^ , ..■,xl, ..., x^"*" ~ , x^"*" ) 
or, equivalently xl_^_i = ^BL{x\), where 0(...) = ^RL is the function of the site values which 
specifies the rule. Here r and / denote the number of cells to the right and to the left of 
site i respectively. The range p = / + r + 1 of the rule is the total number of sites in the 
neighborhood. There are p'' = p^+'+i different configurations on the table for each possible 
combination of site values. Each configuratioun in the table is then simply given by the integer 
number n which runs between and p'^+'+i — 1 (then tolerance = 1/2). They compare to the 
dynamical configuration reached by site i and its r and / first-neighbors at time t and given by 
Ylk=-rP''~^^^)^ ■ "^^^ latter is the input of the rule. The outputs a^ for each configuration n 
are also integers E [0,p — 1]. An integer number R can then be given in base 10 to fully specify 
the rule Rp as R = z2n=o o,nP^- With all these correspondences the following expression is 
obtained 

x^i = Yl «n^ n - ^ /+'■xj+^ - (3) 

n=0 \ k=—r J 

Eq. ([3]) describes all first-order-in-time deterministic CA rules in ID with no freely adjustable 
parameters: the p^+'+i coefficients a„ directly specify the dynamical rule. For example, for Wol- 
fram's rule ""^1102 (ao,ai, ...a-j) = (0, 1, 1, 1,0, 1, 1,0) (see Fig. S2 in the Supporting Information 
[23], where all above notations are also clarified). 

Since each site on the ring satisfies Eq. ([3]), the whole set of equations is globally invariant 
upon translation modulo N^. As shown in the Supporting Information [23], Eq. ([3]) is also 
invariant under certain permutations of the integers in [0,p — 1] (change of colors), reflection, 
change of base, change of range, and shift (Galilean invariance) (see Figs. S3-S5 in the Support- 
ing Information [23J). These invariances allow to classify all CA rules into equivalence classes 
reducing enormously the number of rules to a fewer representative ones. Polynomial maps in- 



volving integers can also be directly derived from Eq. ([3]). For example, the map for all 256 
Wolfram rules ^R^ is 

xi+^ = ao(l - xi+i)(l - xl){l - x\-') + arx\-\l - 4+'){l - x\) + a^x^l - x\+'){l - x^') 

+aexi-^'xi{l - xl~') + ajxl^'4xl' (4) 

where all variables and parameters involved are either '0' or '1' (since p = 2), see Theorem Sll 
in the Supporting Information [23j. Thus, rule ^1102, has the remarkably simple expression 
a^t+i = xl + x\~ {1 — x\ — x\x]^ ). Many Wolfram rules can be analytically solved for the 
orbit so that the solution at each site value i is known in closed form as a function of the 
initial condition and time. These orbits are the CA analogy for the trajectories in differential 
calculus. For example, rule ^122 has orbit x\ = Xq — x^^x^^ so that each site value i is known 
in closed form as a function of the initial condition and time. Then, since rule ^342 is related 
to rule ^122 through a shift transformation, the orbit of the former is automatically known as 
x\ = Xq~ — Xq" Xq" '^ . All 4 trivial local rules '^i?2 (see Fig. S6 in the Supporting Information 
[23]) and all 32 rules ^R^ and ^B^ (which implement the 16 logical functions, see Figs. S7 to 
S9 in the Supporting Information |23j ) can be solved for the orbits as well. Tables SI to S4 
in the Supporting Information [23j show the polynomial maps for all these simple rules and 
the orbits which solve the maps for all rules '^-R2! ^^-^2 ^^'^ ^-^2- Remarkably, rule "62 (and its 
class equivalents ^62, "^92, ^92) is, although predictable, far more complex than the rest of the 
p = 2 rules (see Figs. S8 and S9 in the Supporting Information [23]). While the latter are 
monotonic (changing the input value of one site makes the output to increase (decrease) or to 
remain constant for every possible input value of the other site), rule "^62 is non- monotonic. It 
implements the addition modulo 2 of the input site values. I call this kind of rules Pascal rules: 
they reproduce, for given p and p the Pascal simplex modulo p (see Figs. SIO and Sll in the 
Supporting Information [23j) and are therefore invariant upon addition modulo p. Since the 
outputs of these rules depend only on the sum of the previous site values, these rules constitute 
also a particular case of the so-called totalistic rules [U [2] (which are themselves a subset of all 
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the rules described by the universal map, Eq. ([3|)). Pascal rules are also a subset of the so-called 
additive CA rules. In |24j an algebraic theory is provided for additive CA. 

Eq. ([3]) allows to prove a theorem (see the Supporting Information [23], Theorem S6) which 
shows how any rule is constructed from rules of lower range. It implies, for example, that two 
rules with p = 2 and vectors (ao,ai, ...,a2r+i+i_i) and {ao,ai, ...,a2r+i+i_i,aQ,ai, ...,a2r+i+i_i) 
describe identical behavior (although the latter has a neighborhood with one site to the left 
more than the former). Thus rules ^62 ^^"^ ^10^2 with vectors (0, 1, 1, 0) and (0, 1, 1, 0, 0, 1, 1, 0) 
respectively, describe the same dynamics. The additional degree of freedom of rule ^1022 can 
be used to break the symmetry of the Pascal rule ^^62 from which it is entirely constructed. If 
the symmetry is broken so that some few more (but not all) configurations within the Pascal 
simplex change their output in a monotonic manner, the most complex rules are generated. This 
is the origin of complexity in ID CA. Class 4 rules can be derived through this process. In the 
case of Pascal rule (0, 1, 1,0), the symmetry can be weakly broken by changing the output of 
configuration 'Oil' from zero to one. This amounts to change 03 —t- 1 in rule ^1022 (see Fig. 

2 and Fig S12 in the Supporting Information [23]) obtaining (0,1,1,1,0,1,1,0) i.e. Wolfram's 
celebrated Class 4 rule ^1102, which is known to be capable of universal computation [T] IT0ll25] . 

Figures 3 (center), and 4, show some other examples of Class 4 rules constructed in this way: 
glider guns (coherent structures that radiate other coherent structures periodically) and other 
elements necessary for universal computation are observed in all cases. All these Class 4 rules are 
derived from the Pascal rule ^270301 with vector (0, 1, 1, 0, 1, 0, 0, 1, 1,0, 0, 1, 0, 1, 1, 0) (see Fig. 

3 left). The triangle of expansion shows that certain configurations yield their output within 
the triangle of expansion and can be tuned breaking non- monotonic turns. Examples are con- 
figurations '1010', '0011', '1111', '0010', '0001'. The rule is then copied to the higher range p = 5 
as (0, 1, 1, 0, 1, 0, 0, 1, 1, 0, 0, 1, 0, 1, 1, 0, 0, 1, 1, 0, 1, 0, 0, 1, 1, 0, 0, 1, 0, 1, 1, 0) which corresponds to 
rule ^17714651102, a Pascal rule, as well, which coincides with the original Pascal rule ^27030| 
as a result of Theorem S6 in the Supporting Information [23]. Now, if any of the above configu- 
rations within the triangle of expansion is considered with the added degree of freedom one has 
the configurations: '01010', '00011', '11111', '00010', '00001', which correspond to numbers '10', 



00001 00001 00001 

00011 00011 00011 

00121 00101 00111 

01331 01111 01101 

14641 10001 11111 

Pascal triangle Pascal triangle Wolfram's rule 

modulo 2 110 

FIG. 2: Scheme showing the relationship of the Class 4 Wolfram's ^1102 rule to the Pascal 
triangle modulo 2. A weak symmetry breaking of the addition modulo 2 within the triangle is 
introduced: All the outputs but one are still dictated by the same prescription that generates 
the Pascal triangle. The only exception is configuration 'Oil' that returns '0' as output in the 
Pascal rule while it returns '1' in Wolfram's rule. This slight difference (indicated in green) 
propagates within the triangle of expansion of the rule and is enough to make the resulting 
pattern both complex and unpredictable. 

'3', '31', '2','1' in the decimal system. Then any of the following positions (or several of them) in 
the vector characterizing the Pascal rule ^17714651102 can be tuned to yield Class 4 behavior: 
aio, as, asi, 02, and ai. Figure S13 shows some of these possible symmetry breaking tunings 
and the resulting Class 4 behavior. This result is remarkable: finding these rules by brute force 
simulations might imply the evaluation of thousands of millions of rules. And with the simple 
prescription given above, these rules can be directly found! If addition modulo p is broken in a 
way that the borders of the Pascal simplex are affected (or that non-monotonic turns are created 
within) randomness and Class 3 behavior occurs instead (see Figs. 3 (right) and Fig. S14 in the 
Supporting Information [23j). 
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FIG. 3: Spatiotemporal evolution, obtained from Eq. ([3|), of Pascal rule ^27030| (left) and of 
the Class 4 rule 21771466134^ (center) and the Class 3 rule ^17714692061 (right) derived from 
the former by a weak and a strong symmetry breaking of the addition modulo 2, respectively. 
A detail of the triangle of expansion for each rule is shown on the top. In the case of the Class 
4 rule, the outputs in the initial time steps are the same as in the Pascal rule and only three 
site values (indicated with light blue, denoting '1') are different at t = 3. In the Class 3 rule 
the symmetry breaking is strong, and 12 sites (indicated with light blue, denoting '1' and dark 
blue, denoting '0') have already different values compared to the Pascal rule at t = 3. 



Another consequence of the Theorem S6 in the Supporting Information |23) is that if two 
rules are monotonic with the same trend, the constructed rule is also monotonic. This allows 
to define a complexity index k for every CA, with values k = 1, 2, 3 depending whether the rule 
is monotonic, non- monotonic, or weakly non-monotonic (weakly breaking the addition modulo 
p symmetry), n = 1 describes Class 1 and 2 behaviors, and k = 2 and 3 Class 3 and 4 



9 





FIG. 4: Spatiotemporal evolution of Class 4 rules obtained from a weak symmetry breaking 
of the Pascal rule ^270301 with vector (0,1,1,0,1,0,0,1,1,0,0,1,0,1,1,0), whose triangle of 
expansion is shown above each rule. In red are indicated the configurations that are tuned 
after the rule is copied to higher range. The rules have following codes: (a) ^17714661342, (b) 
^17714661421, (c) ^1771466136^ and (d) 23918949782^. 
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behaviors, respectively. This complexity index is tabulated for all Wolfram rules in Table S5 
in [23]. This complexity index should not be confused with the one introduced by Chua in 
|14j . which is related to the complexity of implementing CA rules by means of cellular neural 
networks. The complexity index introduced in this preprint is directly related to the complexity 
of the dynamical behavior observed and can be calculated for each rule and Theorem S6 in the 
Supporting Information [23j by decomposing the rule in its constructing sub-rules layer by layer. 
The complexity index applies not only to boolean rules with p = 2, r = 1 and / = 1, but to rules 
of any range and any number of symbols in the alphabet. 

Eq. ([3]) can be generalized to higher dimensions and topologies (see Fig. S15 and the 
accompanying discussion in the Supporting Information p3j). For example, the popular Con- 
way's "Game of Life" gS E7| is given by the map xH^ = B (3 - Yll,m=-i4'^'''^^"'^l) + 
x^'-'i3 ( 4 — ^^ m=-i -^t '' ™' 2 ) '^here x]'-' denotes the central site of a Moore neighborhood in 
2D. When p — )• 00 any local map involving real numbers can also be approximated by ^Rp rules 
as briefly shown in the Supporting Information [23J for the case of the logistic equation. 

The ideas and methods presented here can be fruitfully used in manifold ways to model 
complex systems and might be helpful in finding accurate coarse-grained models of physical 
systems with many degrees of freedom [28] . The role played here by the Pascal rules is analogous 
to the one played by separatrices between nonlinear resonances in Hamiltonian systems. The 
destruction of KAM tori in Hamiltonian systems through weak perturbations, creating chaotic 
layers that coexist with islands of regular motion [2^, seems to be mimicked, at the CA level, by 
the weak symmetry breaking of the addition modulo p reported here leading to Class 4 behavior. 
Research on these intriguing connections might also yield valuable insight in the understanding 
of turbulence. 
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SUPPLEMENTARY INFORMATION 



OVERVIEW OF THE ARGUMENT 



A mathematical function can be thought as a rule such that, given a numerical input, the 
output is read from a table containing a series of configurations to which the input compares. 
The complete table gives the rule, i.e. what to do with an appropriate input. Usually, the 
configurations on a table must be sharply and unambiguously distinguished from each other 
within a certain tolerance. Then the table provides the output corresponding to the configuration 
that successfully matches the input. It is then said that the output is the result of the rule acting 
on the input. 

The above paragraph can be seen as a computational way of thinking Mathematics. In the 
following, the argument is inverted. Mathematics is taken as starting point and I shall consider 
the most elementary mathematical structure that imitates computation as described above. Each 
computational process is considered here as an "experiment" and I use standard mathematics to 
theoretically model such an experiment. Proceeding in this way is advantageous since a general 
mathematical theory for cellular automata (CA) can then be rigorously established. 

I consider first the most simple rule, i.e one containing just only one configuration with one 
nonzero output. When input coincides with configuration within a certain tolerance, the rule 
returns output. Otherwise, the result is zero. I seek now for the simplest mathematical structure 
where all actors involved in the above computational process enter. It can be formally depicted 
in the following way: 

output X B {con figuration — input, tolerance) (5) 

The function B has two arguments and, in order for the latter expression to return output 
or zero, it is clear that B should return 1 or 0, respectively. Specifically, it should return 1 
if {configuration — input\ < tolerance and otherwise. It is clear then that B is an even 
function of its first argument. The mathematical function B{x,e) having all these properties is 
the boxcar function plotted in Fig. SI and introduced in Definition SI below. 
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B{x,^) 



Figure SI: The boxcar function B{x, e) is equal to one for — e < x < e and zero otherwise. 



Definition SI: The boxcar function B{x,e) is defined, for e and x real, as 



B{x,e) 



1/x + e x — e 



2 \\x + e\ \x — e\ 



where \x\ 



X X > 
—X X < 



is the absolute value function. □ 



(6) 



The boxcar function is defined over the real numbers and has a discontinuity at x = =t|e|. 
When e = 1/2, the boxcar function is usually called rectangular function. In Appendix A several 
alternative definitions of B and some useful results -some well known- are outlined. 

When one has a rule with several different configurations indexed by n so that to 
configuratiorin corresponds an outputn, the rule is completely given by summing over all ele- 
mentary structures describing the action of the rule on each separate configuratiorin as given 
in Eq.Q. Each configuration is exclusive, i.e., the input cannot simultaneously be equal to two 
different configuratiorin within tolerance. The output of the rule is then 

output = y^ outputn X B{configurationn — input, tolerance) (7) 

ndtable 

The tolerance is related to the distance separating two adjacent configurations on the table and 



is defined as 



tolerance 



configurationn+i — configurationn 



(8) 
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and so, if each configurationn is given by an integer number, tolerance = 1/2 (if each con- 
figuration is instead given by a rational number separated a distance 1/d from the next, then 
tolerance = l/{2d)). Remarkably, and as shown in Section II, a CA can be always made to fit 
in the general structure given by Eq. ([7]). 

II. GENERAL THEORY FOR CELLULAR AUTOMATA (CA) 

A. Universal CA map 

A spatially extended system consisting of a lattice of discrete sites is now considered. The 
topology and dimensions of this lattice can be arbitrary. I focus first on a ID ring containing a 
total number of Ng sites, and develop the complete theory for such a situation. 

An input is given as initial condition in the form of a vector xq = (xq, ...,Xq *). Each of the 
Xq is an integer in the range through p — 1. The subindex labels the discrete time step in the 
evolution and the superindex specifies the position of the site on a ID ring. At each time step t 
the vector xj = (xj, ...,Xj ^) specifies the state of the CA. Inputs and outputs from the rule are 
always integers on the range through p — 1. Periodic boundary conditions so that x^ "^ = x\ 
and x^ = Xf.\ Let xj_,_^ be taken to denote the value of site i at time step t + 1. Formally, its 
dependence on the values at the previous time step is given through the mapping 

^t+l — 'r\^t i^t , ■■■iXf-, ...,Xf ,Xf. ) (9) 

where </>(...) is the function of the site values which specifies the rule. Here r and / denote the 
number of cells to the right and to the left of site i respectively. The range p of the rule is 
defined as the total number of sites involved in the rule and is therefore given by p = / + r + 1. 
Since each site i can have p different values (labelled as integers between and p — 1) there is a 
total of rj = p^ different configurations on the table for each possible combination of symbols. 
For a range p there exists a total number oiT = p = pp^ different rules that can be defined in 
this way. Each configurationn in the table is given by integer number n which runs between 
and Vt — \ (and hence tolerance = 1/2) and which is given as a function of the site values in the 
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table X* as 

I 

Y^ /+ra;i+fc (10) 



n = 

k=—r 



The X* here have also values between and p — 1 and specify all possible configurations. The 
above n specifies each static configuration given on the table. They compare to the dynamical 
configuration nt reached by site i and its r and / first-neighbors at time t 



nt 



E /+'^i+' (11) 



The latter is the input of the rule. 



The outputs ttn for each configuration n, given by Eq.(lO), have integer values between and 
p — 1- An integer number R can then be given to specify the code of the rule in the following 
way 

n=0 

i.e. the a„ are the coefficients that accompany the powers of p when writing R in base p. Since 
the rule depends not only on the a„ but also on the number of symbols p involved and on 
the number of neighbors to the left / and to the right r, to define an arbitrary CA rule in an 
unambiguous way I introduce the quantity 

Rp 

This latter quantity will symbolically denote in all the following the same as (f) in Eq. ([9]) above. 
The labels /, r and p accompanying R specify the rule completely and Eq. S can be written 
in a compact and unambiguous way as 

xj+i = 'i?;(xD (13) 

with the understanding that the rule depends not only on x\ but also on the first r neighbors to 
the right and the first / neighbors to the left. For example, Wolfram's rule 110, involves p = 2 
symbols, and one left and one right neighbors, i.e. r = I = 1 (see Fig. S2). Such rule is then 
written as ^1102. 
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Figure S2: Wolfram's rule ^1102: (top) The table of the rule with each configuratiorin (n) and 

each outputn (an) and code R; (bottom) spatiotemporal behavior of the rule running from a 

single initial seed on a ID ring (time runs from top to bottom and space is plotted 

horizontally, under periodic boundary conditions). 



Since each site i G [1,-^s] on the ring satisfies separately Eq. (13), the whole set of A'^, 



equations is globally invariant upon translation modulo Ng, i.e. upon making the change i — )• i+h 
to every site i (where h is an integer number). 

The following set of correspondences between the quantities here defined and the ones 
introduced in Eq. ([T]) in Section I hold 
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configuratiorin 
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outputn 

output 



nt= ^ p 



k+r i+k 




Sr~^ ^k+r^i+k 



X 



t+l 



and the universal map governing the dynamics of ID cellular automata, by using Eq. (JTl), has 
then the form 




(14) 



or, equivalently, by using Definition SI 



n-i 



X 



t+i 




1 I v^' 

2+^-L.k-. 



\ + r^-T!k- 



r,k+T i+k 1 

•P -^t 2 



^ + Ei 



,p 



,k+r^i+k 



•P -^t I 



l-r^ + El 



-r P -^t I 



(15) 



Although Eq. ( 14 ) seems rather complex, it is to be noted that it describes all first-order- 



in-time deterministic CA maps in ID. The map contains no freely adjustable parameters: the 
an directly specify the dynamical rule. When most an are zero, the sum reduces to just a few 
terms. Eq. (14) is the normal form for all T = p = pp^ deterministic ID CA. 



Sometimes, the configuration is specified not in terms of its binary code but over a sum 
carried over all sites in the neighborhood. The output is made in such cases dependent only 
on the sum of the previous values. These CA are called totalistic and they are a subset of the 



total possibilities described by Eq. (14). Each configuration^ in the table is labelled in such a 



case with an integer number s which runs between and G = p{p — 1) (the latter value is the 
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maxinium value that the sum over the site values can attain) 

I 
E ^^^' (16) 



I 
s = 

k= 



The X* have as before values between and p — 1 and specify all possible configurations of the 
symbols. The above s specifies each static configuration given on the table. They compare to 
the dynamical configuration st reached by site i and its r and / first-neighbors at time t 



I 

st=Yl ^^' (17) 



k=—r 

which is the input of the totalistic rule. 



The outputs (Ts for each configuration s, given by Eq.(16), have integer values between and 
p — 1 like the inputs and the output of the rule. An integer number R can then be given to 

specify the code of the totalistic rule in the following way 


R = Y,asP' (18) 

A totalistic CA rule is therefore labelled as 

^rt; 

where the label T is added to avoid confusion with the coding in the normal form. The following 
set of correspondences between these quantities and the ones introduced in Eq. ([7| in Section I 
holds now 



input <^=^ St = /, 






k=—r 

configurationn <^=^ s = > x*^ 

\ k=—r / 

tolerance <^=^ 1/2 

OUtputn <^=^ CTg 



output 



4+1 



And therefore, the universal map governing the dynamics of ID totalistic cellular automata, is 
obtained by using Eq. ^ as 
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^t+1 


s=0 


sB 
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->: 

k=-r 


^i+k 


1 



(19) 



Eq. (19) is a particular case of Eq.(14). Starting from a totalistic rule with vector 



((To, cri...o"e) described by Eq. (19) the vector specifying the normal rule as described by Eq. 



(14) (oo, ai...af2_i) can be calculated from the following expression 



"Ei = -r.P'' + ''^'+'' ~ ^Ei = -r. a;' + * 



(20) 



The following important theorem establishes that the universal CA maps, Eqs. (14) and 



(19) are invariant upon a bijective application that sends the set of p integers {0, 1, ...,p — 1} to 
any other permutation of the same p integers, provided that inputs, output and specifications 
in the table are subjected to the application. An analogy is provided by computer simulations, 
where to visualize the spatiotemporal pattern of the CA a color to each integer number is 
associated. The invariance under change of colors warrants that any palette can be chosen if 
the distinctness of the symbols is preserved, the latter having no effect on the dynamics (the 
same pattern is obtained albeit with different colors). This, of course, must be the case since 
the dynamics must be the same regardless of the labels that used for the symbols. 



Theorem SI (Invariance under change of colors): The universal CA map, Eq. (I4) 
remains invariant after the following set of transformations 



J+k 



4+1 



J+k 



4+1 



+ Y,X^Bi x\ 



m=0 
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m. 



n ^ n = n 



m=0 ^ ^ 

k=—r m,=0 

^"^ / 1\ 

hn' = a„/ + ^ XmB (on' - m, -\ 



(21) 
(22) 
(23) 

(24) 



where the Xm 's are integer numbers (not necessarily positive) so that the transformations between 
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the sets of p integers {0,l...,p — 1} implied by Eqs. (21), (22) and (24) are bijective so that 
I p-i 



k=—rin=0 

when both sides are different to zero 



m,-\-B[x\ 



i+k 



m 



+ E ^'^" \^t 



Vk _ ^i+k 



(25) 



k=—r 



Proof: By inserting Eqs. (21) to (24) in Eq. (14) 
p-i 



^t+l + /_^ ^mB 
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On-m, 
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Bin- Yp 



k+r i+k 
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k=—r 



where to get to the penultimate equahty I have used that the apphcation between the two sets 



of integers {0, l-..,p — 1} is bijective satisfying Eq. (25) and the sum has been reordered. By 
using result (vii) from Appendix A and the fact that xl_^_^ = a„ for the n = J2k=~rP''~^^^''~^^ 
that satisfies n = Ei=-r p''^'' 4'^'' ^ 



Y^rnB 
m=0 



and, therefore, 
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(26) 



(27) 



n=0 



k=—r 



proving thus the invariance of Eq. ( 14 ) upon change of colors. The proof for Eq. ( 19 ) is similar. 



D 



When one considers specific rules, described by the universal maps above, they usually 
break this symmetry as shown in the following example. 

Example: The asymmetric rule ^4583 (see Fig. S3 - top left) with three symbols '0','1' 
and '2' has range p = r + / + l = 2 (and hence Cl = p^ = 3'^ = 9) and coefficients 
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^3837? 



"'lid -3 ..1"^M 

"■'■i.-d:^ .d 



°15605j 





Figure S3: Spatiotemporal evolution of rules '^4583, "^156053, ^SSS?!] and ^170693- Rules on the 

top are related to the ones at the bottom through a change of colors. Rules on the left are 

related to the ones on the right through reflection. 

(ao, ai, ...as) = (0,1,2,1,2,1,0,0,0) (i.e. the number 453 in base 3 in inverse order) in Eq. 



(14). The neighborhood contains only the cell that is updated in the next time step and the first 



neighboring site to the right. I derive now the rule that operates with the symbols 2, 1, in the 
same way as rule ^^4530 does with symbols 0, 1, 2 respectively. Clearly, Aq = 2, Ai = 0, A2 = —2, 



and from Eqs. (23) and (24), the transformed rule {bo,bi, ...bs) = (2,2,2,1,0,1,0,1,2) is ob- 



tained, which corresponds to rule "^ 156051 (see Fig. S3 - bottom left). The colors corresponding 
to symbols with values (blue) and 2 (red) are exchanged in both figures, while sites with value 
1 (green) remain unchanged. 

Rules "^4533 and '^156053 belong then to the same equivalence class under change of colors. 
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If p is the number of symbols in a rule, there are at most p\ rules in its same equivalence class 
(note that there are rules that may be invariant upon change of some colors: in general, the 
cardinal x of ^'^ equivalence class is, therefore, 1 < x ^ pO- 



Theorem S2 (Invariance under reflection): The universal CA map, Eq. {I4) remains 
invariant after the following set of transformations 

(28) 

(29) 

(30) 

Proof: This theorem can be checked directly by making the corresponding transformations in 
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Eq. (14). n 



Invariance under reflection and under change of colors can be followed after the other in 
either direction since both commute. The dynamical behavior of all these rules is closely related 
and all of them belong to the same equivalence class, which, when added all rules obtained 
under change of colors and its reflections, contains no more as 2pl elements. In Fig. S3 it is 
shown how rules "^4533 and ''156O53 are related through reflection to rules ^38373 and ^170693 
respectively. All these rules belong, therefore, to the same equivalence class. 



Theorem S3 (Invariance under change of base): The universal CA map, Eq. (I4) remains 
invariant after the following transformations 



p ^ p {p > p) (31) 

I 
n ^ n' 

k=—r 



I'^Yl P"^'^'^' «n ^ ^n' = an (32) 



Proof: It can be checked directly by making the transformations implied by Eqs.(31) to (32) in 
Eq. ([14]). D 
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Example: Rules ^1102 and ^5906013 are equivalent after a change of base and exhibit the 
same dynamical behavior. 



The validity of the following theorem can be again checked directly employing Eq. ( 14 ) 



Theorem S4 (Invariance under change of range): The universal CA map, Eq. (I4) 
remains invariant after the following transformations 



I -^ I' {I' > I) 

(r' > r) 
n' = n + a 
br,i = a„ \/a 



r ^ r 



n 



'"n ' "n' — "n 



where 



-r-\ V 

k=-r' k=l+l 



i+k 



(33) 
(34) 
(35) 
(36) 

(37) 



Example: Rule "^62 has vector (00,01,02,03) = (0,1,1,0). If now the range of the rule is 
increased adding one site from the left, i.e. Z — t- /' = / + 1, then il' = p +''+^ = Qp = 8 and a. 



from Eq. (37), can be either or 4. Therefore, the class equivalent rule after change of range is 
(ao, ai, 02, 03, 00, ai, 02, 03) = (0, 1, 1, 0, 0, 1, 1, 0) i.e. rule -'^1022. 



Example: Rule ^62 has also vector (oq, oi, 02, 03) = (0, 1, 1, 0). If now the range of the rule is 
increased adding one site from the right, i.e. r — )• r' = r + 1, then il' = p'+'' +^ = Qp = 8 and 



(T, from Eq. (37), can be either or 1. The class equivalent rule after change of range is then 
(00,00,01,01,02,02,03,03) = (0,0,1, 1, 1, 1,0,0) i.e. rule ^QQ\. 



Theorem S5 (Invariance under shift): The universal CA map, Eq. (I4) remains invariant 
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after the following sets of transformations 
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and 



l-l 



r + 1 



X 



J-i 

t+i 



I -^ l + l 
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(38) 
(39) 
(40) 

(41) 
(42) 
(43) 



Proof: By making the transformations implied by Eqs.(38) to (40) in Eq. (14) 



j;{ = x;a.sfn- g /'+'^+ixr',j') = i:««^f--E^' 



+r i—l-\-k 

Jb-t- 



1 



(44) 

?i=0 \ fe'=— r— 1 / n=0 \ k=—r / 

where the change of the dummy variable A;' — )• /c — 1 has been made. The latter expression is, of 



course, equivalent to 



'^t+i ~ ^p[^t ) 



(45) 



Making the change z — >• i + 1 (since global translation invariance holds) the invariance under 



shift is proved. The proof involving Eqs. (41) to (43) is similar. □ 



The invariance under shift of the universal CA map implies the existence of classes of rules 
related by a breaking of this symmetry which contain, at most, p elements. The rules on these 
classes share the same code R but the neighborhoods contain different numbers of sites to the left 
and to the right (although the range p is the same). These rules satisfy the following identities 

(46) 



Rp\Xi] 



^-'r;+H4-') 



'^'i^'iA"'] 



(47) 



The rules are equivalent in a sense that, when known the dynamical behavior of one of them, 
the behavior of the others is predictable in terms of the latter just by applying a global spatial 
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shift of the dynamical state on the ring. 

Example: Rules "'^1102, ^1102 and "^llOl (see Fig. S4) are related by a shift transformation 
and belong to the same equivalence class under shift: once known the dynamical behavior of 
rule ^1102, the ones of rules ^1102 and '^llO^ are equivalent by globally shifting the dynamical 
state on the ring one site to the right or or to the left, respectively. 

Rules related by shifting transformations correspond to the same dynamical behaviors as 
seen by observers moving with different (constant) velocities to either side of the ring. For a 
given rule 'i?L the shifted rules '•^'"Klp^'" corresponds to the same dynamics as followed by an 
observer moving on the ring at a constant velocity if (where the positive sign corresponds to 
motion to the left). This is the principle of Galilean invariance for cellular automata rules. 

Shift transformations and invariance under change of range are particularly powerful 
when combined together to relate symmetrical rules (i.e. with / = r) with the same range 
p but a different value of R. For example, rules ^42 and ^42 are trivially related by a shift 
transformation. Since rules ^^42 and ^682 are equivalent under change of range, as are the rules 
^42 and ^482, it is clear that rules "'^482 and ^683 are also related by a shift transformation. This 
is not trivial at all, since ^482 and ^682 are both symmetrical rules, with the same range, but 
with a different value of R. Both of them represent, however, the same dynamics as followed 
by observers moving at different relative velocities. This is clearly evidenced in Fig. S5 where 
the spatiotemporal evolution of all these rules is shown. 

Theorem S6 (Constructor's theorem): Let a set of p rules each denoted by '(^m)^ with 

m G [0,p — 1] with code Am = S„=o ^n p"", with Q = |)'+'"+i. A rule '+^i?p can be constructed 

from the left in the following way 

p-i 

'2 



'+'r; =Y,B ( x^+'+i - m, i ) \Am)l (48) 



m=0 
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^110° 




Figure S4: Spatiotemporal evolution of rules ^1102, ^1102 and ^^llOl. All three rules are related 
by a shift transformation and belong to the same equivalence class under shift. Rules related 
by shifting transformations correspond to the same dynamical behaviors as seen by observers 
moving with different (constant) velocities to either side of the ring. For a given rule Rp, the 
shifted rules RL^"" corresponds to the same dynamics as followed by an observer moving on 
the ring at a constant velocity zizv (where the positive sign corresponds to motion to the left). 
This is the principle of Galilean invariance for cellular automata rules. 
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Figure S5: Spatiotemporal evolution of the rules indicated in the figure and their 
interrelationships through shift transformation and change of range. 
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Proof: By using Eqs. Eq. p^ and (|48) 
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has been introduced. By using now {xii) from Appendix A 
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where n' = X]fc=-rP '"^ 

rule, Eq. Jsil has then the code R = Y.n>=l Cn'P^' = En=o Emio ai™V"+'^"- 

The proof of the statement on the construction from the right proceeds in a similar manner. 



By using Eqs. Eq. (14) and (49) 
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has been introduced. By using now (xiii) from Appendix A 
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'+^^+2 = pJ7 and Cn' = a)"7 ,, . The 



new rule, Eq. (51) has then the code R = Yln'=o ^n'P"' = X]n=o X]m=o '^" ^"^^™ ^^'^ ^^^ 
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result is proved. □ 

Example: Theorem S6 can be straightforwardly applied to any rule. Wolfram Rule ^302 
is known to be a random number generator. It has vector (00,01,02,03,04,05,06,07) = 
(0, 1, 1, 1, 1,0, 0, 0). To obtain the construction from the left, simply separate (0, 1, 1, 1, 1, 0, 0, 0) 
into p = 2 consecutive parts with same size. Rules with vectors (0,1,1,1) and (1,0,0,0) are 
obtained, which correspond, respectively, to rules ^^142 and "^Ig- To construct the same rule 
from the right, separate the odd entries and the even entries of (0,1,1,1,1,0,0,0). Rules 
(0, 1, 1, 0) and (1, 1, 0, 0) which correspond to rules ^62 and ^82 are hence obtained. □ 

Definition S2 (Copying rules): The rule ^AL is said to be copied to a higher range if it is 
used as the only rule in constructing '^^i?2 from the left, i.e. if\Am)p = '^2 ^"^ ^'^ ^Q- {4^)- ^ 

Example: To copy rule ^62 to a higher range p = 3 its vector (00,01,02,03) = (0,1,1,0) 
is simply concatenated to itself, adding the same sequence of zeroes and ones after it i.e. 
(0, 1, 1, 0, 0, 1, 1, 0). The latter corresponds to rule ^1022 which is the copy of "62 to range p = 3. 
This is consistent with the construction from the left of rule ^1022 as established in Theorem S6. 
Both rules belong to the same class under change of range and display an identical dynamical 
behavior. D 

Example: To copy rule '^263 to a higher range /) = 2 its vector (00,01,02) = (1,2,2) 
is concatenated after itself two times (since p = 3) i.e. (1,2,2,1,2,2,1,2,2). The latter 
corresponds to rule ^189253 which is the copy of the original ^^253 rule. □ 

The following theorem allows to construct a polynomial of integer variables for each CA rule 



contained in Eq. (14) and Eq. (19) and provides, indeed, the universal CA maps in polynomial 
form. This is an important result since it allows to pass from a description in terms of boxcar 
functions to polynomials involving integer variables. In order to obtain the latter a linear 
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system of equations involving a Vandermonde matrix of 17 x ri dimensions needs to be solved. 



Eq. (14), which is computationally inexpensive and sufficient to characterize each CA rule is, 



therefore, more fundamental. 



Rli^l) = A+\ = X] "" Z] v'^ ' ■ ^t 



Theorem S7 (Polynomial maps): The following result holds 
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where the coefficients an 's (rational numbers) are the solutions of the following linear system of 
equations 
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(53) 



'rt;{x\) = X, 
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(54) 



where the coefficients a^ 's (rational numbers) are the solutions of the following linear system of 

equations 
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(55) 



Proof: Eqs. (52) and (54) are, respectively, the polynomial interpolations of Eqs. (14) and 



(19). For the topic of polynomial interpolation see for example [30] □ 



[30] M.J.D. Powell (1981). Approximation Theory and Methods, Chapter 4. Cambridge University Press. 
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B. Results for boolean {p = 2) CA rules 

The simplest CA rules involve two symbols p = 2 and only the cell xl that is updated by 
the rule (r = I = 0). Since the input can be just only zero or one, and for each posibility the 
output can be also only zero or one there is a total of four such rules, labelled as ^R2- 

Theorem S8 (Universal CA map for the 4 most simple boolean rules): The universal 



CA map, Eq. (14), reduces to the following expression 



xi^^ = ao{l - xi) + aixi (56) 

for the 4 most simple boolean rules (p=2, 1=0, r=0), i.e. the rules ^R2, where R = X]„=o^"^"- 



Proof: By inserting / = r = 0, p = 2, fi = 2 in Eq. (14) 
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aoB{xl,-j +aiB(xi-l,- 



whence, by simply applying (xv) from Appendix A, the result is proved. An alternative proof 



is given by using Theorem S7, Eq. (52), which takes in this case the simple form 



x^+i = ao + oiXt 



(57) 



where ao and ai can be obtained by solving Eq. 53 



1 

1 1 






= 


ao 
_ai_ 



(58) 



which yields ao = ao and ai = ai — uq thus proving again Eq. (56). □ 



The maps for the ''i?2 rules are summarized in Table SI and plotted for an initial condition 
consisting on a single seed (Fig. S6, top) and a random initial vector (Fig. S6, bottom). 

The local and global dynamics of the repetitive application of these rules is trivially pre- 
dictable: given a vector xo = (xj, ..., Xq ^) as initial condition, the vector x^ = (x|,...,Xj ") can 
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Figure S6: Spatiotemporal evolution of the four rules '^i?2 (values of R are indicated over each 
panel) for an initial condition consisting on a single seed (top) and for a random vector 

(bottom). 
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Table SI: The four cellular automata maps for the simplest boolean rules i?2 
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Table S2: Local (top) and global (bottom) dependence in time of the state of the cellular 
automata for the simplest boolean rules '^i?2 

be determined as a function of xq and t by means of a closed analytical expression. In Table 
S2, these expressions are given for the local (top) and the global (bottom) dynamics. 

One can easily prove all entries on Table S2. The proof of the entry for rule ^^12 proceeds by 
induction. Clearly, the expression is valid at t = and at t = 1. If one assumes it valid for t — 1, 
then, at time t 



1 



H-l 



1 



(-1) 



'-'^h + lil 



-1 



it-1^ 



-(-l)*-Vo + ^(l + (-l)*-^) 



(-l)*4 + -(l-(-l)0 



q.e.d. 
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i-1 
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Table S3: The sixteen maps for the ceUular automata implementations of the sixteen logical 

functions. The maps correspond to rules ^R2- 



Theorem S9 (Universal CA map for the 16 boolean logical functions): The universal 



CA map, Eq. (I4), reduces to the following expression (boolean logical functions) 



J+i 



ao{l - x\){l - x\ ^) + aix\ \1 - x\) + a2x\{l - x\ ^) + a2.x\x 



J^i-l 



(59) 



for all 16 rules with p=2, 1=0, r=l, i.e. the rules R2, and to 



J+i 



ao(l - xr')(l - 4) + aixlil - xl+') + a2xj+^(l - x]) + a^x'+'x 



i+l\ 



J+'i{ 



J+l^i 



for all 16 rules with p=2, 1=1, r=0, i.e. the rules i?2- In both cases R = Yln=o'^n'^"- 



(60) 
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Map 
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4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 
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.*+i 



+ X 



Xf 






i+1 



.i+1 



XlX 



i+1 



t-^t 



Xt 



x\ + xj+^ - 



2x1x1 



i+l 



XI 



^' 



J+^^i 



X 



.i+l 



i+l 



+ xlx 



i+l 



t-^t 



rytl _|_ rr"^ /y^l ^ 



xi + X 



i+l 
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Table S3 (cont): The sixteen maps for the cellular automata implementations of the sixteen 
logical functions. The maps correspond to rules ^i?2- 



Proof: By inserting / = 0, p = 2, r = 1, O = 4 in Eq. (|14|), 

3 

^t+1 



J2anB(2xi + xi-'-n,^\ 

n=0 ^ ^ 



and by using the definition of n given by Eq. ( 10 ), i.e. n = 2x^ + x* ^, the latter expression can 

be rewritten as 

11 

^~l^ 1 



1 1 

a;'=Oa;»-i=0 



iB{2{xl-x') + {x 



i-1 



"2 
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By applying now result (xii) from Appendix A 

1 1 



x»=Ox»-l=0 



«-^ ^. ._-,.,. -.„,„,., -.-,1 



The sum can now be carried explicitly by replacing the corresponding values of x* and x"^ ^ which 
pick values either or 1 yielding four terms. Then, by applying result (xv) from Appendix A, 



Eq. (59) is proved. 



The proof of Eq. (60) is similar after making the replacements xj — )■ xj^ , xj 



X* — 7- X*'*'"'", X* "*"—>• X*. □ 

The rules ^R2 and ^R2 which implement each of the 16 logical boolean functions are shown 
in Fig. S7 together with their associated Venn diagrams. For both types of rules '^i?2 ^'^d ^-^2' 
A denotes the input on site x] (which is updated at the next time step) while B denotes the 
input on site xl~ for rules ^R^ or x^"*" for rules ^i?2- 

As an example the logical function B NOT A, returns 'TRUE' (i.e. 1) if and only if in the 
previous time step only the site B (and not A) has the value 'TRUE'. This means for rule "^22 
that one has the truth table xj^^ = 1 (i.e. 'TRUE') only if A is 'FALSE' (i.e. if xj = 0) and B 
is 'TRUE' (i.e. x^~ = 1). Rule ^42 implements the same logical function, but now B 'TRUE' 
means x]^ = 1. 

Because these functions shall prove important in building the theory of complexity presented 
here, I give in Table S3 the full analytical expressions obtained from Theorem S9. In Figs. S8 
and S9, the spatiotemporal evolutions of these rules are plotted. 

Although some of the patterns seem rather complex, their behavior is entirely predictable, 
as established in the following theorem that I give without proof (all results can be proved by 
induction -by using the results in Appendix A- and can be easily checked with a computer). 
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Figure S7: The sixteen boolean logical functions implemented through the cellular automata 

rules ^R2 and ^i?2 indicated in the figure. The associated Venn diagrams are also shown. For 

both types of rules '^i?2 ^^^ ^-^2' ^ denotes the input on site xl (which is updated at the next 

time step) while B denotes the input on site xj^ for rules ^R^ or xj for rules ^i?2- 

Theorem SIO (Predictability of all boolean rules ^R^ and ^R2)' If^o = {x^, ...Xg, ...Xq ") 
denotes the initial state of a CA at time t = 0, for time t > 1 (positive integer), all rules ^R^ 
and ^i?2 ^''"^ predictable and the value at the site on a later time t, xl, can be known for each 
rule, as a function oft and the initial site values xq. The orbits (local dependence in time for 
the state of the cellular automata) for all rules ^R^ and ^R2 are given in Table S4- 
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Figure S8: Spatiotemporal evolution of the sixteen rules ^R2 {R is indicated over each panel) 
for an initial condition consisting on a single seed (top) and a random initial vector (bottom). 
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Figure S9: Spatiotemporal evolution of the sixteen rules ^i?2 (^ is indicated over each panel) 
for an initial condition consisting on a single seed (top) and a random initial vector (bottom). 
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Table S4 (cont.): The orbits (local dependence in time for the state of the cellular automata) 

for all rules ^R^ and ^R^. 

Specially interesting are the rules "^62, ''92, ^62, ^92- Although predictable, the behavior of 
these rules is far more complex than the one of the other rules. They correspond to the logical 
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XOR function and to the complementary XNOR function. In the foUowing section, it is shown 
how this is the building block of complexity in physical systems. These rules represent addition 
modulo 2 of both site values at a previous instant of time. Therefore, they are also totalistic 
rules. Rules °6^, ^6^, °9^ and ^9^ have totalistic codes °2^2^ ""IT^, ^5T^, ^5T^ respectively. 

Theorem Sll (Universal CA map for the 256 boolean Wolfram Rules): The universal 



CA map, Eq. (14), reduces to the foUowing expression 

x^+i = ao(l-xj+i)(l-x^)(l-xr)+aixr(l-:rj+i)(l-xj) + 

+ a2xlil - x^+i)(l - xj-i) + asxlxl'(l - x\^^) + a,x\^\l - x\){l - xj"') + 

+ a5xj+^xj-i(l - x\) + a6xj+^xj(l - x\-^) + a^x'^^^x\x\-^ (61) 

for all 256 rules with p=2, 1=1, r=l (Wolfram rules), i.e. the rules ^R^, where R = X]n=o'^"^"- 



Proof: By inserting / = r = 1, p = 2, fi = 8 in Eq. (14) 



t:j+i = Yl ^nB (4xj+i + 2xi + xr 1 - n, ^] 



(62) 

n=0 



and by using the definition of n given by Eq. (10), i.e. n = 4:X^~^^ + 2x* + x* ^, the latter 



expression can be rewritten as 
1 1 1 



X 



t+i 



E E E «4.>+i+2..»+x.«-i3(4(xr-^'+') + 2(xj-x^) + (xr^-x-i),^] 



a;*+l=0a;»=0x»-i=0 

By applying now result (xii) from Appendix A 
1 1 1 

Ei+l=0a''=0x»-i=0 

The sum can now be carried explicitly by replacing the corresponding values of x*"^^, x' and 
x*~^ which pick values either or 1 yielding eight terms in the sum. Then, by applying result 



E E E «4x.«+^+2.»+.'-i^(^^'-^*+\^)^(^^-x\^)e(xri-x-\^] 



(xv) from Appendix A, Eq. (59) is proved. □ 



The 256 specific maps obtained from Theorem Sll are provided in Table S6 in Appendix B 
for each rule (the values for all a„ specifying each rule are there in each case explicitly indicated) . 
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Theorem S12 (Some equivalence classes of boolean rules): Let '^2 ^'^'^ '^2 ^^ ^^^ 
different boolean rules, with A = l^„=o '^"^" '^^'^ ^ ~ Sn=o ^«2"' ™^^ ^ ~ 2'+*^'+^. and 
J2k=-r 2^~^^x^^^ The following results hold 



n 



• {i) B'2 is related to vlg through a change of colors transformation (also called global 
complementation) when bn = ^ — an-i-n for all n. 

• {a) ^ B2 is related to A2 through reflection when bn = fl„_y-' (2k+r_2^->')x^+'' f^^ ^^^ 
n (e.g. for all x^^^ , where k G [—r,l]). 

• {Hi) ^-62 ^-5 related to A^ through reflection following a change of colors transformation 
when bn = I- a^^.^.^+j-! __ (2fc+r_2i-fc)2;i+fe for all n (e.g. for all x'^^ , where k e [-r, I]). 

Proof: This theorem is indeed a corollary of Theorems SI and S2, for p = 2. The only bijec- 
tive application to change the colors that can be defined is a; — )■ 1 — x which coincides with GC. □ 



Example: The above theorem can be cursorily checked with the specific case of Wolfram 
^1102 rule for which fi — 1 = 7, (00,01,02,03,04,05,06,07) = (0,1,1,1,0,1,1,0). Therefore, 
global complementation (from (i) in the above theorem) gives (60, ^1, &2, ^3, &4, &5, ^6, ^7) = 
(1 — 07, 1 — 06,1 — 05, 1 — 04,1 — 03, 1 — 02,1 — ai,l — oo) = (1,0,0,1,0,0,0,1), which 
corresponds to rule ^1372- Under reflection one obtains, from {ii) (60, ^1, &2, ^3, &4, ^5, &6, &?) = 
(00,04,02,00,01,05,03,07) = (0,0,1,1,1,1,1,0) which corresponds to rule ^1242- Fi- 
nally, global complementation after reflection gives, from (Hi): (60,^1,^2,63,^4,65,66,67) = 
(1 — 07, 1 — 03, 1 — 05, 1 — oi, 1 — 06, 1 — 02, 1 — 04, 1 — Oo) = (1, 0, 0, 0, 0, 0, 1, 1) which corresponds 
to rule H93^. n 

Example: Consider now ^63, for which. $7 — 1 = 3, (oq, oi, 02, 03) = (0, 1, 1, 0). Global comple- 
mentation (from (i) in the above theorem) gives (60, 61, 62, 63) = (1 — 03, 1 — 02, 1 — oi, 1 — oq) = 
(1,0,0,1), which corresponds to rule '^92- For the rule related through reflection one obtains, 
from (ii) (60,61,62,63) = (00,02,01,03) = (0,1,1,0) which corresponds to ^62 (note that / 
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and r are exchanged). Finally, global complementation following reflection gives, from (in): 
(60, ^1, b2, b^) = (1 — 03, 1 — oi, 1 — a2, 1 — ao) = (0, 1, 1, 0) which corresponds to rule ^92 (note 
again that / and r are exchanged compared to the original rule ^62). □ 

Theorem S12 allows the 256 Wolfram rules to be classified in 88 equivalence classes. The 
shift transformation combined with change of range (see Theorems S4 and S5 and Figure S5) 
allows to reduce this number even further. Finally, carefully applied and combined with the 
constructor's theorem S6, Theorem SIO (Table S4) allows even analytical expressions for the 
spatiotemporal evolution of some of the rules to be easily found, since they are also Wolfram 
rules after being copied to range p = 3. Analogous expressions as in Theorem SIO (Table S4) 
might be found for most of the Wolfram rules and this problem will be addressed elsewhere. 

Rules "'^1102 and ^542 as well as their class equivalent rules, display a more complex behavior 
than the rest of the Wolfram rules (of range p = 3). It was observed above that rule "62 and 
its class equivalent rules display also a more complex behavior than any other boolean rule of 
range p = 2. Rules ^1102 is even more complex than rule ^62 because it is unpredictable. In 
the following section a general theory of complexity is introduced, allowing the wide variety of 
behaviors found in CA to be identified. 

C. General theory of complexity 

Based on extensive computer simulations, Stephen Wolfram introduced a classification of CA 
dynamics in terms of observed behaviors of increasing complexity [31] Although this classification 
was empirical, based on computer experiments, its robustness is compelling. As shown below, 
the theory of CA presented here elucidates from a fundamental (mathematical) point of view 
this classification, and provides a rigorous definition of complexity for the dynamics of cellular 



[31] See Wolfram, S. Cellular automata as models of complexity. Nature, 311, 419 (1984) and Wolfram, S. 
Universality and Complexity in Cellular Automata. Physica D, 10, 1 (1984). 
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automata. The theory presented here is vahd for rules of any range and involving an arbitrary 
number of symbols. These ideas might be extended to continuum systems: in the last section I 



show how maps involving real numbers arise from the universal CA map (Eq. 14). 

Speaking about a generic initial condition, Wolfram classified CA dynamics into four main, 
broad classes, as a function of the limiting behavior reached by the CA 

• Class 1: Evolution leads to a spatially homogeneous state 

• Class 2: Evolution leads to a set of separated simple stable or periodic structures 

• Class 3: Evolution leads to a chaotic, aperiodic or nested pattern 

• Class 4: Evolution leads to complex, localized structures, some times long-lived 

This classification, although being entirely descriptive and qualitative, is consistent with a 
huge number of computer simulations. Some of the rules exhibit a certain coexistence between 
classes. 

I introduce now some definitions and theorems to proceed further. The proof of Theorem 
S13 is trivial from the constructor's theorem S6. 

Definition S3 (Monotonic rules): A local rule ^Rp is monotonic of sign + / - if increasing 
the value of xl makes the value of xl_^i to increase/ decrease. If it remains constant, the rule is 
called neutrally monotonic. A general rule ^RL is monotonic, if it is a monotonic function of 
each of every site values in the neighborhood. 

Theorem S13 (Monotonicity from construction): Let a set of p rules each denoted 
by \Am)p with m G [0,p — 1]. Then ^^^KL, constructed from the left, and ''KL'^^, con- 
structed from the right, are monotonic if all constructing rules {Am)p are monotonic of the 
same sign for every site value (some of the constructing rules can also be neutrally monotonic). 
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Definition S4 (Non-monotonic turns in construction layers): A rule RL is said to 
have a non-monotonic turn in the construction layer k if in the construction process starting 
recursively from the monotonic rules ^lil, the rule *" R'^ with l'-\-r' = k<l + r, constructed 
after k steps, is non-monotonic. The higher the number of non-monotonic turns, the stronger 
the non-monotonicity of the rule. 

Every rule ^i?2 ^-^d ^R2, with the exception of "^62 and its class equivalents, is monotonic. 
This can be seen from the constructor theorem. Rules ''i?2 ^-^^ constructed from rules '^i?2 adding 
one site to the right. Because all rules '^-R2 ^^^ monotonic, any construction involving rules with 
the same monotonicity sign is monotonic as well, as a consequence of Theorem S13. The only 
exception is the construction involving simultaneously rules "^12 and "^22: since the former is 
monotonic of sign - and the latter monotonic of sign +, the constructed rules are not monotonic. 
The latter rules construct together ^62 and its class equivalents! These rules have, therefore, a 
non-monotonic turn on the first layer. One can easily see that rule "^62 is nan — monotonic: let 
the site on the right have value '0'; then, increasing the site on the left from '0' to '1' makes the 
output also to increase from '0' to '1'; however, if the site on the right has value '1', increasing 
the site on the left from '0' to '1' makes the output to decrease. 

The analysis of the most simple boolean rules in Theorems S9 and SIO and Figures S6, S8 and 
S9 showed indeed that rule ''62 and its three class equivalents upon reflection and change of colors 
is more complex as any other rule ''i?2- And what the analysis indeed reveals as well is that, 
since additional degrees of freedom cannot destroy the non-monotonic character achieved by rule 
"62, Class 1 and Class 2 cellular automata correspond to monotonic rules, and the most simpler 
non-monotonic rules begin with the nested structure, which corresponds already to Class 3. This 
elucidates an observation made by Wolfram [32] that the number of Class 3 rules increase with 
increasing range: for an increasing number of construction layers, there are an increased number 
of possibilities to create non-monotonic turns. And when a non-monotonic turn is created, any 



[32] See Wolfram, S. A New Kind of Science, Wolfram Media (2002), p. 948. 
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rule constructed after that by employing the previous one, is also automatically non-monotonic 
because non-monotonic turns in previous layers cannot be destroyed. 

Another important feature of such a rule as ^62 is that it puts into play every available symbol 
starting from a single seed. Such rules, yielding nested structures, are our starting point in our 
aim to understand the emergence of complexity since, being regular and predictable, they lie at 
the verge of complexity and allow to classify every other rule upon perturbations and symmetry 
breaking. I show this below. 



I first note the following important invariance satisfied by the totalistic CA map, Eq. (19) 



and state it as a theorem. This invariance can be directly checked on Eq. (19). 



Theorem S14 (Invariance upon addition modulo p): The totalistic universal CA map, 



Eq. (19) remains invariant after the following set of transformations 

s — >• s' = s + mp (63) 

Us — )■ O"^/ = (Ts Vs (64) 

where m is an integer number so that s + mp £ [0, p{p — 1)]. 

Most rules break this symmetry. The only exceptions are, of course, the ones that satisfy 
o"s = (Ts+mp, V?7i so that s + mp G [0, p(p — 1)]. And specially interesting within these are 
those which put every symbol [0,p — 1] into play during the time evolution. Because of their 
importance, and for reference, I name these rules invariant upon addition modulo p Pascal 
rules, after the French mathematician Blaise Pascal: these rules, as observed below, reproduce 
after a bijective application all Pascal simplices modulo p. When p = 2, the Pascal simplex 
coincides with the Pascal triangle. When p > 2 the Pascal simplex is related to the multinomial 
expansion modulo p. Pascal rules are a subset of the so-called additive cellular automata, for 
which an algebraic theory was formulated. [33] For reference, I introduce these rules as a 

[33] See Martin, O.; Odlyzko, A. M. and Wolfram S. Algebraic properties of cellular automata. Comm. 
Math. Phys. 93, 219-258 (1984). 
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definition. 

Definition S5 (Pascal rules): A totalistic rule '■RTp is called Pascal rule if it satisfies the 
following property 

(^s+mp = CTs = S (65) 

where m is an integer number so that s + mp S [0, p{p — 1)]. These rules perform the addition 
modulo p of all site values contained in the neighborhood. 

Theorem S15 (Pascal rules perform the addition modulo-p): A Pascal rule RTp per- 
forms the addition modulo p of all site values contained in the neighborhood, i.e. it has the 
form. 

m=0 s=0 \ k=—r / \k=—r / 

where TZp{x) is the remainder upon division of an integer number x by p. 



Proof: The first of the equahties in Eq.(66) is obtained from the totalistic universal map Eq. 



(19) by making the transformation s — )■ s-\-mp and using the definition of the Pascal rule above. 
Then, since = p{p — 1) 

x\^, = Y^a^BL- ^xj+^M 

s=0 V k=-r / 

p(p-l) / I 



Y^ (js+mpB is+mp- Y a^t■^^ 2 ) 

s+mp=0 \ k=—r / 

p p-l / I \ 



m=0 s=0 \ k=—r / 

The second equality is then proved by using result (xvi) from Appendix A, with Up = p.O 

The structure of these rules is pretty simple. They have vectors (cJo, ...,o"pp) with the 
structure (5,5,5"), where S is a chain of integers 012. ..{p — 1) repeated until the pp positions 
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• 



• 



characterizing the rule are fiUed. S" is the chain S truncated when position p{p — 1) is reached. 

Examples: 

02ri with (do, ai, as) = (0, 1, 0), i.e. xj+i = 7^2 (eL-i 4'^'' 
HOT^ with (ao,ai, (72,^3) = (0,1,0,1), i.e. xj+i = 7^2 (eL-i^'^^ 
HOTi with (ao, (71,(72,^3,^4) = (0,1,0,1,0), i.e. x^^^ = 7^2 (EL-i^^t^'' 

0l02r3l with (C70, (71, (72,(73, (74) = (0,1,2,0,1), i.e. 4+i = 7^3 (EL-1 ^^t^'' 
^588r3^ with ((70,(71,(72,(73,(74,(75,(76) = (0,1,2,0,1,2,0), 

i.e. x|+, = 7^3(EL-l^^')• 
°1346680^5^ with ((70,(71,(72,(73,(74,(75,(76,(77,(78) = (0,1,2,3,4,0,1,2,3), i.e. i.e. xl_^_^ 



I546268555T5I with ((70, (71,..., (711, (712) = (0,1,2,3,4,0,1,2,3,4,0,1,2), i.e. xj+i 
°275781750T6i with ((7o,(7i, ...,(79,(710) = (0,1,2,3,4,5,0,1,2,3,4), i.e. xj+i 
°78050441406r7i with ((7o, (7i, ..., (7ii, (712) = (0,1,2,3,4,5,6,0,1,2,3,4,5), i.e. xj+i 



All Pascal rules fall in the Wolfram class 3, since they yield nested structures that can be very 
complex (see Fig. SIO), but that are regular and with predictable features: since they reproduce 
the Pascal simplices, they are related analytically to the multinomial expansion. The structures 
that are formed from a single seed are contained in a big triangular structure. The latter is 
called triangle of expansion of the Pascal rule. For rules with p = 2, the triangle of expansion 
coincides with the Pascal triangle modulo p. Figure Sll shows a detail of the first time steps of 
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^lOT^ 




^588TJ 




''27578175076 








"13466801' 




°78050441406T' 



i n a A a 6 a 



i J t a ^ ^ ^ J a ^ J a ^ 



Figure SIO: Spatiotemporal evolution of some Pascal rules with totalistic codes indicated in 
the figure. The number of different colors in the figures coincides with p. Regular, nested 

structures arise in every case. 
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Figure Sll: Detail of the former stages of the spatiotemporal evolution from a single seed of 

the Pascal rule with totalistic code ''l346680r5 . The Pascal structure modulo 5 is clearly 

recognized. The '1' colored red at the bottom of the figure corresponds to the value 6 in the 

Pascal triangle, which happens to be equal to 1 modulo 5. 

the evolution of rule ''l346680r5 showing how, indeed it calculates the Pascal triangle modulo 
5. 

Rule ^62 is also a Pascal rule with totalistic code ''2T2^ (see Fig. SIO). Any other rule ^R^ 
(with exception of the class equivalent rule ^Og) is simpler than the Pascal rule: either the 
triangle of expansion disappears, or it is homogeneously filled, as it is the case with rule '^142- 
This is the consequence of the absence of the non-monotonic turn, characteristic of rule '^62. 

It is, indeed, possible, to have more complex structures than the Pascal rules. When a 
Pascal rule is copied to a higher range and then only a few configurations within the triangle of 
expansion are tuned so the non-monotonic turn of the copy is destroyed, the symmetry of the 
rule is broken introducing defects that cause coherent structures to arise and propagate always 
within the triangle of expansion. This is, as shown below, a weak symmetry breaking of the 
addition modulo p. Any other way of breaking the symmetry of the Pascal rule is stronger and 
merely reduces or conserves the complexity of the reference Pascal rule. I state now the central 
result of this whole research and substantiate it with examples. 



Class 4 CA arise from a weak symmetry breaking of the addition modulo p 
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Example: Starting with the Pascal rule ^2T2, it is first converted to the normal code, which, as 
known from above is "^62 (in general, Eq. (20), should be used). The rule has vector (0, 1, 1, 0). 



Its spatiotemporal evolution starting from a simple seed shows that the only configuration of 
two sites contained within the triangle of expansion yielding its output within the triangle of 
expansion is '11', i.e. '3' in base 2. Rule ^^62 is then copied to the higher range /) = 3, obtaining 
(0,1,1,0,0,1,1,0), i.e. rule ^1022 and the output of configuration 'Oil' falling within the triangle 
of expansion is switched. This amounts to merely make the change 03 — )• 1 in rule ^1022. In 
this way the addition modulo p symmetry of rule *'2r2 is weakly broken and, at the same time, 
the non-monotonic turn of the first substructure of the rule is removed, i.e. (0,1,1,0,...) (non- 
monotonic) — )• (0,1,1,1,...) (monotonic). Then (0,1,1,1,0,1,1,0), is obtained, i.e. rule ^1102 which 
is the celebrated Wolfram rule, known to be a universal Class 4 CA. Figure S12 sketches this 
process. From the point of view of the construction process, ^1102 is non-monotonic because it 
contains rule "63 in its second construction layer. But is not so strongly monotonic as rule ^1022 
because the latter is entirely constructed by rule "62 and has, therefore two non-monotonic turns. 
Rule ^1102 contains also the monotonic rule '^142 in its second construction layer. Rule ^1102 
is weakly non-monotonic in contrast with the strongly non-monotonic Pascal rule '^62 (or ^102^). 

Example: Consider the Pascal rule ^10r|. First, the rule is converted to the normal 



code by using Eq. Q and rule ^27030^ with vector (0,1,1,0,1,0,0,1,1,0,0,1,0,1,1,0) is 
obtained. The triangle of expansion shows that certain configurations yield their output within 
the triangle and can be tuned breaking non-monotonic turns. Examples are configurations 
'1010', '0011', '1111', '0010', '0001'. The rule is then copied to the higher range p = 5, 

(0,1,1,0,1,0,0,1,1,0,0,1,0,1,1,0,0,1,1,0,1,0,0,1,1,0,0,1,0,1,1,0) 

which corresponds to rule ^17714651102. Now, if any of the above configurations within the 
triangle of expansion is considered with the added degree of freedom one has the configurations: 
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Figure S12: Scheme showing how a weak symmetry breaking of the Pascal rule *'2r2 (= ^62) 

leads to Wolfram Class 4 rule ^1102- (a) Triangle of expansion for rule ^2T2 at the former 

stages of growth, (b) In range p = 2, only the configuration '11' belongs to the triangle of 

expansion yielding its output '0', within the triangle of expansion, (c) The rule is copied to 

range p = 3 yielding the equivalent rule ^1022 with an additional degree of freedom that can be 

used to break the Pascal rule on the configuration '11' (d). The resulting rule ^1102 has an 

unpredictable pattern within the triangle of expansion as a consequence of the weak symmetry 

breaking of rule '^2r2^ introduced by switching 03 from '0' to '1' in rule ^1022- 

'01010', '00011', '11111', '00010', '00001', which correspond to numbers '10', '3', '31', '2','1' in the 
decimal system. Then any of the following positions (or several of them) can be tuned to yield 
Class 4 behavior: oio, 03, 031, 02, and ai. Figure 4 in the main text shows some of these possible 
symmetry breaking tunings and the resulting Class 4 behavior. This result is remarkable: 
finding these rules by brute force simulations might imply the evaluation of thousands of 
millions of rules. And with the simple prescription given above, these rules can be directly found! 

Example: Also totalistic rules can weakly break the addition modulo p. In fact, the following 
example illustrates more precisely the difference between weak and strong symmetry breaking. 
While the former eliminates a non-monotonic turn when breaking the modulo p symmetry, the 
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Figure S13: Spatiotemporal evolution of Class 4 rules obtained from a weak symmetry breaking 

of the totalistic Pascal rule ^WT^ {=^ 27030^ with vector (0,1,1,0,1,0,0,1,1,0,0,1,0,1,1,0)), 

whose triangle of expansion is shown above each rule. In red are indicated the configurations 

that are tuned after the rule is copied to higher range. The rules have following codes: (a) 



2l771466134|, (b) 2l771466142|, (c) ^17714661361 and (d) 23918949782 



2- 



56 



latter does not. As a consequence the Pascal triangle structure is severely distorted and a 
mixing of all the available symbols ends in chaotic behavior. This latter behavior is less complex 
than the one obtained through the, more subtle, weak symmetry breaking. It is nonetheless 
important, that such complex rules can be engineered using the Pascal simplices modulo p as 
template. 

Figure S14 illustrates the difference between weak and strong symmetry breaking of a Pascal 
rule. The Pascal rule has vector (o"o,o-i,(T2,o"3,(T4,(T5, ue) = (0,1,2,0,1,2,0) meaning that the 
output value is each time a situation in which the sum of the site values is either 3 or 6. These 
positions are responsible for the breaking of the monotonicity of the rule since the natural drive 
of the sum to increase in the natural numbers is reseted to zero at the sum values which are 
multiples of 3. The symmetry of the Pascal rule can be now broken by switching the value 
of one of the inner configurations in the triangle of expansion. If the innermost one, o"g, is 
switched it is clear that a structure similar to the Pascal rule propagating unchanged during 
several time steps is formed. Then a defect occurs and the behavior of this defect decides wether 
the rule behaves then "chaotically" or with a Class 4 behavior, yielding complex structures. 
There are two possibilities to destroy the symmetry of the Pascal rule: by tuning a^ to '1' or 
to '2'. The first tuning keeps the non-monotonic turn of the Pascal rule while at the same 
time destroying the modulo p addition symmetry. The rule is therefore strongly non-monotonic 
and the resulting rule, with vector (cjo, (Ji, (T2, (T3, CJ4, CJ5, ug) = (0, 1, 2, 0, 1, 2, 1) yields seemingly 
chaotic Class 3 behavior. If now the innermost configuration is tuned as ag — ?• 2, the modulo p 
addition is destroyed together with the hard non-monotonic turn, and the gentle monotonicity 
of the previous steps saturates to a value of two. This illustrates the weak symmetry breaking 
of the addition modulo p. The resulting rule (o"o,cri,o"2,o"3,o"4, (T5, de) = (0,1,2,0,1,2,2), has 
Class 4 behavior. In fact, this seems characteristic of all observed totalistic rules of Class 4: the 
saturation of the sum within the triangle of expansion meaning that, after a certain transient 
(that can be very long) no complex structures survive anymore. 

It is to be noted that the whole discussion about complexity in this section applies to all 
members in a given equivalence class under change of colors, reflection, shift, etc. 
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Figure S14: Spatiotemporal evolution of Pascal rule with code ^SSSTg^ (top). The modulo p 

symmetry can be broken so that the strong non-monotonicity of the Pascal rule is preserved. 

In such case rule ^ISlTTg (bottom left) is obtained and the corresponding dynamical behavior 

is chaotic (Class 3). If however, the symmetry breaking is weak, so that the modulo p 

symmetry is destroyed together with the non-monotonic turn, the Class 4 totalistic rule 

^2046T3 (bottom right) is obtained. The vectors of each totalistic rule are indicated on the 

figure, and in red are marked the positions that break monotonicity within the Pascal simplex. 
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Rule 
iQi 

I2I 
131 

141 

1 

2 

1 

2 



I5I 



i6i 



171 
i8i 



191 
ilQi 
illi 

1121 

1 
2 



I13I 



I14I 
I15I 



ii«i 



18;^ 



I10I 



19i 



lool 



^22, 



I23I 



I24I 



(00,01,02,03 
(0,0,0,0 
(1,0,0,0 
(0,1,0,0 
(1,1,0,0 
(0,0,1,0 
(1,0,1,0 
(0,1,1,0 
(1,1,1,0 
(0,0,0,1 
(1,0,0,1 
(0,1,0,1 
(1,1,0,1 
(0,0,1,1 
(1,0,1,1 
(0,1,1,1 

(1,1,1,1 

(0,1,0,0 
(1,1,0,0 
(0,1,1,0 
(1,1,1,0 
(0,0,0,1 



04,05, 


26,a7j GG 


LK 


UULK 


K 


0,0,0,0; 


12551 


loi 


1255^ 





0,0,0,0; 


11271 


'4 


il27^ 




0,0,0,0; 


11911 


II6I 


1247^ 




0,0,0,0; 


I63I 


1171 


III9I 




0,0,0,0; 


1223^ 


141 


1223^ 




0,0,0,0; 


I95I 


I5I 
^2 


I95I 




0,0,0,0; 


1159^ 


I20I 


1215^ 




0,0,0,0; 


I31I 


121° 


i87i 




0,0,0,0; 


1239^ 


i64i 


1253^ 




0,0,0,0; 


lull 


I65I 


I125I 




0,0,0,0; 


11751 


i80i 


1245^ 




0,0,0,0; 


1471 


i8li 


I1I7I 




0,0,0,0; 


1207^ 


i68i 


1221^ 




0,0,0,0; 


I79I 


I69I 


I93I 




0,0,0,0; 


I143I 


i84i 


1213^ 




0,0,0,0; 


I15I 


I85I 


I85I 




1,0,0,0; 


II83I 


II8I 


1183^ 




1,0,0,0; 


I55I 


I19I 


I55I 




1,0,0,0; 


11511 


I22I 


I151I 




1,0,0,0; 


I23I 


I23I 


I23I 




1,0,0,0; 


1231^ 


i66i 


1189^ 


2 



TABLE S5 : The 88 independent rules under global complementation (GC), left-right 

transformation (LR) and global complementation of the left right transformation (GCLR) out 

of the total of 256 rules ^R^ and their class-equivalent rules and complexity index. 
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Rule 



I27I 



loal 
1 

2 

1 

2 

1 

2 



28 



I29I 



I30I 



I32I 



I33I 
I34I 
I35I 



iQAl 



^36 



I37I 



I38I 



140,1 



I41I 
I42I 
I43I 



1/1 /ii 



144 



I45I 



i46i 



I50I 



I51I 
I54I 



(.00,01,02 

:i,o 

0,1 

1,1 
0,0 
1,0 
0,1 
0,0 
1,0 
0,1 

1,1 

0,0 
1,0 
0,1 
0,0 
1,0 
0,1 

1,1 

0,0 
1,0 
0,1 
0,1 

1,1 

0,1 



03 

0,1 





1 



04, 
1,0 
1,0 
1,0 
1,0 
1,0 
1,0 
0,1 
0,1 
0,1 
0,1 
0,1 
0,1 
0,1 
0,1 
0,1 
0,1 
0,1 
0,1 
0,1 
0,1 

1,1 
1,1 
1,1 



05,06,07) 

0,0) 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 



GC 



11031 
11671 



11991 
1711 



11351 

'2511 

11231 
11871 

'591 

i219i 



1911 



11551 

1235^ 

I107I 

11711 

1431 

i203i 



I75I 



I139I 
I179I 



LR 

1671 
i82i 
'831 
I70I 



I71I 



'861 
I32I 



I33I 



IaqI 



1491 



iQfil 
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I51I 
I147I 
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I29I 
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1 

2 
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'2411 
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i89i 



1209^ 
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I51I 
'147I 



K 

2 
2 
1 
2 
1 
2 
1 
2 
1 
1 
2 
2 
2 
2 
2 
1 
1 
2 
2 
1 
1 
1 
2 



TABLE S5 (cont. 
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'154I 
1156^ 
'I6OI 
'1621 
'1641 
'1681 
'iTQi 
'1721 
'1781 
'1841 
1200^ 
1204^ 
i232i 



(.00,01,02 

(0,1 
0,0 
0,1 
0,0 
0,1 
0,1 
0,1 
0,0 
0,1 
0,0 
0,0 
0,1 
0,0 
0,0 
0,1 
0,0 
0,1 
0,0 
0,0 
0,0 
0,0 



03 

1,0 





1 



04,05, 


36,07J LrCv 


LK 


UULK 


K 


o,o,o,r 


'158^ 


'148^ 


'2I4I 


2 


o,o,o,r 


1238^ 


'192^ 


i252| 


1 


o,o,o,r 


I174I 


1208^ 


1244^ 


1 


o,o,o,r 


1206^ 


'196^ 


i220| 


1 


o,o,o,r 


il42i 


'212^ 


'212^ 


1 


i,o,o,r 


'182^ 


il46| 


'182^ 


2 


i,o,o,r 


'150^ 


'150^ 


'150^ 


2 


1,0,0,1^ 


1230^ 


'194I 


'188^ 


2 


1,0,0,1^ 


'166^ 


'210^ 


'180| 


2 


1,0,0,1^ 


'198^ 


'198^ 


il56| 


2 


0,1,0,1^ 


i250| 


il60| 


1250^ 


1 


o,i,o,r 


'186^ 


'1761 


1242^ 


1 


o,i,o,r 


'218^ 


II64I 


'218^ 


2 


o,i,o,r 


1234^ 


1224^ 


1248^ 


1 


0,1,0,1^ 


'170^ 


i240| 


i240| 


2 


0,1,0,1^ 


1202^ 


i228| 


'216^ 


1 


1,1,0,1^ 


'178^ 


'1781 


'1781 


1 


1,1,0,1^ 


i226| 


1226^ 


'1841 


1 


0,0,1, r 


i236| 


1200^ 


1236^ 


1 


0,0,1, r 


1204^ 


i204| 


i204| 


1 


0,1,1,1; 


1232^ 


i232| 


1232^ 


1 



TABLE S5 (cont. 
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I introduce a complexity index k with values 1, 2 or 3, according to the following prescription 

• K = 1: Neutrally monotonic and other monotonic rules (Class 1 and Class 2 behaviors) 

• K = 2: Strongly non-monotonic rules, Pascal rules and non-monotonic rules that strongly 
break the addition modulo p (Class 3 behavior) 

• «; = 3: Weakly non-monotonic rules that weakly break the addition modulo p (Class 4 
behavior) 

This index allows to classify every CA rule. As an example, in Table S5, the complexity 
index k is listed for the 256 Wolfram classes such as they fall in the 88 equivalence classes 
under change of colors and reflection. (As noted above, shift invariance also implies that not 
all of these 88 rules are independent since rules 12 and 34 are also equivalent). The index is 
calculated easily from the construction theorem. 

Example: Wolfram Rule ^SO^ has vector (00501,02,03,04,05506)07) = (0,1,1,1,1,0,0,0). To 
obtain the construction from the left simply separate (0, 1, 1, 1, 1, 0, 0, 0) into p = 2 consecutive 
parts with same size. Rules (0,1,1,1) and (1,0,0,0) are obtained, which correspond, respec- 
tively, to rules "142 ^^^ ^^2- '^^e latter rules are constructed from the right by rules (0, 1) and 
(1, 1) (Rule "^142) and (1,0) and (0,0) (Rule '^12)- The rules in the first construction layer are 
both monotonic, but they have different monotonicity sign: rule '^142 has monotonicity sign -|- 
while rule "I2 has monotonicity sign -. As a result, rule ^302 has a non-monotonic turn in its 
second layer and belongs to either Class 3 or Class 4. Compared to the closest Pascal rule "62 
copied to higher range, i.e. °1022, one sees that 03, 04, 05 and oe are switched. The positions 
controlled by 04, 05 and ag affect the central spine of the triangle of expansion of the associated 
Pascal rule, and break strongly the addition modulo p. Rule ^302 belongs, therefore, to Class 
3 and has complexity index k = 2. □ 
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D. CA in higher dimensions 

The above results can be easily generalized to an arbitrary number of dimensions. In 2D, for 
example, by using the site that is updated after each time step as a pivoting site and defining ID 
codes for each spatial direction counterclockwise, universal maps for all possible deterministic 
CA in 2D (depending on the topology of the interactions in the lattice) can be derived. The 
most popular neighborhoods in 2D are shown in Fig. S15. The von Neumann neighborhood 
can be specified giving two binary codes ni and n2 as indicated in the figure. The hexagonal 
neighborhood requires the specification of three codes to describe each configuration. Finally, 
the Moore neighborhood requires four ID codes. In each neighborhood, the site on the center 
is updated on the next time step. 

Since a configuration is now composed of several ID configurations that must hold simultane- 
ously, a product of boxcar functions is needed to fully account for the dynamical state. For a von 
Neumann neighborhood with horizontal and vertical ranges given by the pairs pi = /i + ri + 1 



and p2 = h + ^2 + I {^i = p^^ and ^2 = p'^^) respectively, it is obtained, instead of Eq. (14) 



^ii = E' E'an.n.^ U - E Z+'-^-i^'^ 2 ) ^ ( "2 - E /+-r^;'^■+^ ^ ) (68) 

ni=0 712=0 y k=—ri j \ ]z=—T2 I 

(note that in the figure only the entirely symmetrical neighborhood with pi = p2 = 2>\s shown) . 

For a totalistic rule over a von Neumann neighborhood 

(pi+P2-l)(p-l) f h h 



E asB{s + l-Y.<^''-Y.-t'^"^\\ (69) 

s=0 \ fc=— n k=—T2 

For an hexagonal neighborhood 

^ii = E E E ««i-->"3^ -1 - E p'^"-^'""', 

ni=0 112=0 713=0 \ k=—ri 

xB (n2 - E p'^'-'xt''''^',\\ bL,- E p'^'''x]:^+', \ I (70) 

\ k=-r2 I \ k=-r3 



and for the totalistic case 

(pi+P2+P3-2){p-l) 






s=0 




E asBis + 2-Y. at''^-' - E 4"'''^' " E -^'^'' ^ (^1) 
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von Neumann 



Hexagonal 



Moore 




(ni,n2)=(1,5) 



(ni,n2,n3) = (3,2,6) (ni,n2,n3,n4)=(4,1,1,5) 



Figure S15: Symmetrical Von Neumann, hexagonal and Moore neighborhoods in two 
dimensions. A configuration can be specified by giving two, three, or four ID codes, 

respectively. 



Finally, for the Moore neighborhood one has, for the general case, 

ni-lQ2-ir23-104-l / h 



4 J 



E E E E «---,n3,n4^ 1-1 - E i^^''^^"''-' 



ni=0 ?i2=0 ^3=0 n4=0 
k=—r2 
k=—r4 



k=—r 



'2 



'2 



k=—rz 



B\n^- Y, p'=+'^^xi+'^'-'+'' 



.. 1 
2 



(72) 
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and for the totalistic one 

(pi+P2+P3+P4-3)(p-l) / li I2 h 

i,j V^ n I I o V^ i+k,j-k Sr^ i+k,j ST^ i+k,j+k 

s=0 \ k=—ri k=—r2 k=—r3 

(73) 



c=-r4 J 



k= 

Above I considered the most general situation, where each spatial direction can be non- 
symmetrical respect to the site i^j. In many interesting cases, however, symmetrical neigh- 
borhoods as the ones in Fig. S15 are considered, when particularized to them, these expres- 
sions reduce to very simple forms. Taking now, for example, a Moore neighborhood with 
li = H = H = Ia = fi = T~2 = T2, = '''A = 1) a-s in Fig. S2b, the latter expression simplifies 
to 

s=0 \ fc,m=— 1 / 

A famous totalistic CA in such a Moore neighborhood due to Gerard Vichniac is the so-called 
"vote" and is given by a vector (ctq, cji, ..., ag) = (0, 0, 0, 0, 1, 1, 1, 1, 1). For such CA, since the 
sum over the cells in the neighborhood can only be at most '9' and the CA only returns '1' when 
s > 5 and zero otherwise, 

s=5 \ k,m=—l j s=5 \ A;,m=— 1 / 

where the sum has been extended to infinity (a trick that is possible to employ here for this 
specific CA because of its structure). By using result (x) from Appendix A 

x% = r{ E xi+'^'^+™-^) (76) 

The r.h.s. of the latter expression is only one when the sum over the cells is higher or equal 
than 5. 

Another famous semitotalistic 2D boolean CA called "Game of life" and invented by J. H. 
Conway is also defined on such a Moore neighborhood in terms of the following rules, governing 
the behavior of the site in the center of the neighborhood (this site is called a "live cell" if it 
has value '1' and a "dead cell" if it has value '0') 
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• 1. Any live cell with fewer than two live neighbors dies, as if caused by under-population. 

• 2. Any live cell with two or three live neighbors lives on to the next generation. 

• 3. Any live cell with more than three live neighbors dies, as if by overcrowding. 

• 4. Any dead cell with exactly three live neighbors becomes a live cell, as if by reproduction. 
This set of rules can indeed be reduced to three 

• 1. If on the entire neighborhood the sum of all site values is '3' the cell in the center lives 
(or becomes a live cell) in the next generation. 

• 2. If on the entire neighborhood the sum of all site values is '4' the cell in the center lives 
in the next generation only if it is already alive. 

• 3. The cell in the center is dead on the next generation if neither 1. nor 2. holds 

By employing Eq. (JTl) an analytical map for this CA can be written directly, translating the 
rules into the mathematical language introduced here 

x^I,=bU- J2 xr"+^M+xi'^-^|4- J2 ^^'■'■'■"^4) (77) 

\ k,m=—l j \ fc,m=— 1 / 

This is a local map for each site i, j in the lattice that implements the Game of Life. 
E. Continuum limit 



Although, as defined, Eqs. (19) and Eq.(14) apply to a finite set of integers G [0,p— 1] which 
are mapped to the integers in the same interval, if p — t- oo the above expression may reproduce 
with arbitrary precision any continuous map. I assume / = r = and a number of symbols p. 
Then Q = pP = p (since p = l + r + l = l). p can now play the role of the precision in the sense 
that a real number is given in terms of multiples of 1/p and then a mapping of the rational 



numbers contained in the interval [0, 1] to itself is performed. This can be seen from Eq. (14) 
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since 



4+1 



Y^ QnB In - xl, -j 



n=0 
p-1 



y^Qn^s 



n=0 



n — XI 



p-1 ' 2(j9 - 1) 



(78) 



where result (iii) from Appendix A has been used. By defining the real quantities n' = n/{p— 1), 
Ul = xl/{p—l) and yl_^^i = xl_^_-^^/{p—l), by taking the limit p — ;■ oo, and by using {v) in Appendix 
A 

lim / dn- — "^^^6 {n - yl) 

p^~ Jo [p -ly 



Vi+i 



(79) 



The latter expression, involving real quantities constitutes the limiting behavior of a CA map 
involving an infinitely large number of integers. The logistic map, for example, can be reproduced 
if 

an'(p-i) = r{p - lfn{l - n) (80) 



since in such case Eq. ( 79 ) becomes 



yU^=ryl{l-y\) 



(81) 



In terms of the original map involving only integers a CA rule ^B^ described by Eq. ( 14 ) with 

n 



dn = r{p - l)n 1 



p — 1 
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and where the code R, defined by Eq.(12) is given by 

f7-i p-1 

R=y^ a„p" = lim V rn(p - 1)(1 - n/(p - l))p" (83) 

^ — ' p— s>oo ^ — ' 

n=0 n=0 

coincides with the logistic map in its dynamical behavior. Simulations show that if the limit is 



not taken and p = 1000, Eq. (78) (obtained from Eq. (14) in this case) already provides an 
excellent approximation to the logistic map. 
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Appendix A: Some results involving the function B{x, e) 

In this Appendix several results involving the function B are outlined. Proofs are skipped 

for most of the statements, although they can be given easily by induction and exhaustion. 

The following definitions of B{x, e) are equivalent to Definition SI 

1 





B{x, e) = - [sign{x + e) — sign{x - 


e)] 


(84) 




= H{x + e)-H{x-e) 




(85) 




= sign{() F ( 1 2 




(86) 




= sign{e) X(_|,|,|,|)(x) 




(87) 


1 


x > 

is the sign function, Hix) = 


1 

< 


x>0 

is the Heaviside 


-1 


x <0 





x<0 



where sign{x) = 

[-1 X <0 

fl xG(-|e|,|e|) 
function, and X{-\e\,\€\){x) = < is the indicator function for the interval 

I otherwise 

i-\el\e\). 

Results: The function B{x, e) satisfies 

(t) Bi-x,e)=Bix,e) 

(ii) B{x, -e) = -B{x, e) 

(Hi) B{ax, e) = B{x, |), for real a. 

(iv) B"'{x, e) = [fj] B{x, \e\), for n non-negative integer. 

(v) lime_^o 2e ~ '^(^)' where 5{x) is the Dirac delta function. 

(vi) J2k=o ^(^ - 2fce> e) = >3{x - {N - l)e, Ne), for < e < 1/2. 

(vii) J2k=o ^(^ - 2^e> e) = 1; for integer x G [-e, (2iV - l)e] and < e < 1/2. 

(via) ^^Jg B [x — k,^) =1 a,ndJ2k=o ^^ (^ ~ ^' s) ~ ^' for integer x £ [—1/2, {2N— 
l)/2]. 
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M Ek=o fik)B{x-k,l)=fix), for integer xe [-1/2, i2N-l)/2]. 
(x) ET=o I3{x - 2ke, e) = H{x + e), for < e < 1/2. 
(x^) Er=-oo >3{x - 2ke, e) = l,forO<e< 1/2. 

(xii) B{an + bm,e) = B{n,e)B{m,e), for any integers a,b,n,m such that \an\ ^ \bm\ 
when either n or m is non-zero and < e < 1/2. (xiib) 7/6 = one then also has 
B{an, e) = B{n, e) for any integer a ^ 0. 

(xiii) B{n, e) + B{m, e) — B{n, e)B{m, e) = B{nm, e), for any integers n, m and < e < 
1/2. 

(xiv) B{n — a, e) + B{n — b,e) = B {{n — a){n — b) , e) , for any integers n, a, b, {a / b) 
andO <e< 1/2. 



• (xv) B{x, ^) = 1 — X and B{x — 1, ^) = x, when x can be only equal to zero or one. 

• (xvi) lZp{a) = Ylrn=Q Sfc=o kB [a — k — mp, j) where TZp{a) is the remainder of the divi- 
sion of a non-negative integer a by p and Up is a suitable upper-bound (positive integer). 

Proof: Results (i) to (iv) follow directly from Definition SI. Result (v) is a direct consequence 



of Eq. (85) and from the standard definition of a derivative of a function of one real variable 



x (with e — )• playing the role of an infinitesimal increment). Result (vi) can be proved by 
induction: For A^ = 1 the result holds. If the result is considered valid for A^ — 1, one has, for N 

N-^l N-2 

^ B{x - 2ke, e) = B{x - 2{N - l)e, e) + ^ B{x - 2ke, e) 

fc=0 A:=0 



B{x - 2{N - l)e, e) + B{x - {N - 2)e, {N - l)e) 



lfx-{2N-3)e x-{2N-l)e x + e x - {2N - 3)e 



2 \\x-{2N - 3)e| \x - {2N - l)e| \x + e| \x - {2N - 3)e| 
1/X + . -(^^-n^\^B(.-(N-l).,m) 



2 \\x + e\ |x-(2A^-l)e| 
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which proves (vi). Results (vii) to (ix) come then from (vi). Result (x) can also be proved 
using (vi) since 



N-l 



y^B{x-2ke,e) = lim V S(x - 2/ce, e) = lim ^(x - (A^ - l)e, A^e) 

k=0 k=0 



1 / x + e x-(2N-l)e\ 1/x + e 
lim - — z = - ■; r + 1 

N^oo2\\x + e\ |x- (2iV- l)e|y 2V|x + e| 



= H{x + e) 

The proof of (xi) proceeds in a similar way. (xii) can be proved by taking into account that both 
the l.h.s and the r.h.s. are one only when both n and m vanish and zero otherwise (since in the 
latter case \an\ 7^ |6m|). The proof of (xiii) proceeds by exhaustion: considering the possibilities 
a) both n and m zero or b) either n or m or both nonzero, the equality is always valid and these 
are the only possibilities. Result (xiv) is a trivial consequence of (xiii). Result (xv) is directly 
obtained from (viii) and can be easily checked from the Definition SI by applying it to both 
possible values of the variable x. 

For a positive integer p, two integers a > k > are said to be congruent modulo p, if 
their difference a — k is an integer multiple of p, i.e. if there exist another integer m such 
that a — k — mp = 0. If such is the case, the r.h.s. of the expression in result [xvi) is equal 
to one only for one integer m in the interval [0, a] called quotient and for one integer A;, 
called remainder. Since the two sums are carried over all possible values for these integers, 
it is clear, by applying {viii) that the r.h.s of [xvi) provides the remainder. Note that 
Up = a would be a rather conservative upper bound for the integers m over which the sum 
is carried and usually a value Up < a can be found, depending of the specific context on 
which the result is applied. A non-negative integer number a is divisible by p when lZp{a) = 0. □ 
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Appendix B: Maps for the 256 Wolfram rules. 
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Rule 



^0^ 



111 



lol 



I3I 



Ul 



lc;l 



^6^ 



I7I 



% 



^91 



ilQi 



I11I 



11^ 



I10I 



12i 



llQl 



13^ 



HAl 



I15I 



H6i 



I17I 



lisl 



18i 



linl 



19^ 



loni 



20^ 



I21I 



(ao,ai,a2 
0,0 
1,0 



0,1 

1,1 
0,0 
1,0 
0,1 

1,1 

0,0 
1,0 

0,1 

1,1 

0,0 
1,0 
0,1 

1,1 

0,0 
1,0 
0,1 

1,1 

0,0 
1,0 



03 

0,0 
0,0 

0,0 
0,0 
1,0 
1,0 
1,0 
1,0 
0,1 
0,1 

0,1 
0,1 

1,1 
1,1 
1,1 
1,1 

0,0 
0,0 
0,0 
0,0 
1,0 
1,0 



04,05,06,07) 

0,0,0,0 

0,0,0,0 



0,0,0,0 
0,0,0,0 
0,0,0,0 
0,0,0,0 
0,0,0,0 
0,0,0,0 
0,0,0,0 
0,0,0,0 

0,0,0,0 
0,0,0,0 
0,0,0,0 
0,0,0,0 
0,0,0,0 
0,0,0,0 
1,0,0,0 
1,0,0,0 
1,0,0,0 
1,0,0,0 
1,0,0,0 
1,0,0,0 



Map 



4+1 



4+1 






J-1 



X 



i+l 



■«-l 



i+l 



+ x'i-xl +x)-xy' + 






t+1 

i 



X 

^t+l 

^t+l 
x" 



J-1 



«+l^'i-l 



X, x 



xlx 



J-1 



t-^t 



2 + 1 

+ Xf XfX 



i-1 



t-^t 



J+l 



t-\- Xf Xf 






XfX 



i™«-l 



J+l 



t-^t 



+ Xf xix 



i-1 



t-^t 



X 



X 



4+1 

i 
t+l 

^J+1 

i 

t+l 

^\+l 



4+1 



4+1 



.«+l 



J+'l^i 



„i-l 



Xf + X 



+ X 

i-1 
t 



i+l^i-1 



J+l^i-1 



2x1x1'^ + 2x'^^x\x 



Irpirri-l 



J+l 



i ?■— 1 

■^t-^t ^t -^t-^t 



r^i^i-l I rJ+l^irJ-1 

XfXf -|- Xf X^Xj 

«+l™J^« — 1 



J + 1 



„i„i—l 



x\ + X 



i+l, 



'i-1 



+ ^'^' 



+ 



i+l ^i ^i—l 



-\--^XfXf Z,Xf X^Xj 



J-l 



i+l^i-1 



X, X 



J+l 



J+l 






i-1 



^«+l„J^i-l 

lib J- JU4-JU 



X 



t+l 



X 



ryi + l^i 
■^t -^t 

J + l 



J-l 



t + l 
^t + l 

x" 

^t+l 
^t+l 
^t+l 



^«+l^i-l 



+ Xf ■ ^Xf 



+ XfX 



j-l 



t-^t 



t-^t 



J+l^i^i-l 



J+l^i 



x'f'^'^Xf + X 



i-1 



J+l^i-l 



rfirfi-l I rf.i+1 rpi rf.i-1 



•^t-^t 



H-^t 



J+l 



t+l 
i 



J+l 



J+l^i 



i+l^i-l 



tJUJ- Ji/f JbM Uj 






i-1 



t-^t 



Xt — X 



i-1 



+ XfX 



i^i—l 



t-^t 



J+l 



r^i+lr^i I rr^-l _ 0™*+-!-^ 



X 



t+l 



4+1 



^f u.f -r u.f 

J+lry.i-1 



i+l^i-l 
t -^t 



x\x\ ^ + 2x 



i+l 



XlX 



i-1 



t-^t 



I ™«+l™i™«-l 



rr.i+1 _i_ ^i _ O^^+l^-i _ rfi+lrpi-l _ rf.i^i-1 _|_ Orri+lrpirri-l 



X 



t+l 



1 _ ^*+l^i 
1 Xf Xf 



^' 



^t ^ -^t 

J+lr 



J-l 



H-^t 



^t-^t 



Table S6: The cellular automata maps for the 256 Wolfram boolean rules ^i?2 ^s obtained from 

Eq. pi) 
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Rule 

I23I 
I24I 



I25I 



i26i 



loYl 



27i 



loQl 



28^ 



looi 



29^ 



iQni 



30^ 



I31I 



I32I 



iQQl 



33i 



iQ^l 



Mh 



1qc;1 



35^ 



^m\ 



I37I 



I38I 



1391 
I41I 



I42I 



(ao,«i,«2 
(0,1 



1,1 
0,0 
1,0 
0,1 

1,1 

0,0 
1,0 
0,1 

1,1 

0,0 
1,0 

0,1 

1,1 

0,0 
1,0 
0,1 

1,1 

0,0 
1,0 



03 

1,0 

1,0 
0,1 
0,1 
0,1 
0,1 

1,1 
1,1 
1,1 

1,1 

0,0 
0,0 

0,0 
0,0 
1,0 
1,0 
1,0 
1,0 
0,1 
0,1 



04,05,06,07) 
1,0,0,0) 

1,0,0,0 
1,0,0,0 
1,0,0,0 
1,0,0,0 
1,0,0,0 
1,0,0,0 
1,0,0,0 
1,0,0,0 

1,0,0,0 
0,1,0,0 
0,1,0,0 

0,1,0,0 
0,1,0,0 
0,1,0,0 
0,1,0,0 
0,1,0,0 
0,1,0,0 
0,1,0,0 
0,1,0,0 



(0,1,0,1,0,1,0,0) 



Map 



.*+l 



j+1. 



+ x\- 2xl^'x\ + xl 



J-1 



2x^^x^^ 



2x\xt^+ 






H-^t 



rf.'i' + '^rpi 



«+l^*-l 



X+ X 









„*+l 



.i+1 






Xj X( -|- X^Xf- 



J-1 



+ 2xlx 



i-l 



7 + 1 

X I -■- /ya *■ /y* 
J- Uj-j-Uj 



j-1 



t-^t 



J+^ 



Xi + X 



i-1 
t 



J+1 i-1 






«+l. 



J-1 



x; — X 



«+i^«-i 



J-i 



X 



t+i 



4+1 



4+1 



4+1 



4+1 



.*+i 



+ x 



2x^+^x^ 



+ XjXt 

«+l^«-l 



Xj X 



t-^t 



-|- Xj XjXj 



^«+i^j 

X/ x + 



J-1 



+ x!x 



i-1 



t-^t 



J+1 



J+l^i 



+ x| - 2x^^"xJ + xj~^ 



2x^+ixri 



J-1 



+ 



+2xj+^xjxr^ 



H-^t 






i+l^«-l 



«+l^*-l 



Xi X 



J+1. 



X^ X^X^ 



j-1 



X. X 



?+l 

/ya " I -*- /ya (■ /ya 



j-1 



X 



t+1 



J+1 



t-^t 



J+1 



+ ~I~ <^+ •^/^ •^ 



*-l I 9^«+l^«-l I ^«^ 



i-1 
t 



9^i+l^i^i-l 
-ZX^ XjXj 



X 



„i-l 



t+1 
^t+1 

x" 



XiX 



J-1 



t-^t 



J+1 



J+1 



J+1. 



/ya t" I /ya "r J- ~, 4 I fya 

Xj -|- X( Xf "I" -^i 



i+l^i-1 



„J+l™i™J— 1 



t+1 
j 



X'' /in ' o^ I /ya 



i+l^i-l 



XiX 



J-1 



t-^t 



J+1 



J+1 



„i-l 



«+l^«-l 



+ 2xr'x: 



j+1 

XI -'- ry* ^ /y> 
X 0-1 + Jb 



j-1 



t-^f 



i-1 



XL /ya ^ /ya '' I /ya 

1 - xj+1 + xj+^xr^ 



J-1 



j+1. 



Oo^'^O^'^ -1 _1_ /ya*' T J- ,-^2 ,ya 



J-1 



t-^t 



XjX 



j-1 



t-^t 



X 



t+1 



4+1 



4+1 



^i+l^i-1 I ™i™J-l 
x^ Xj -|- XjX^ 



9^«+l^i™*-l 
Z,X^ XjX^ 



J+1 



x! + X 



i+1. 



„i-l 



J+l^J-1 



+ 2xJ^'xr + 2xjx: 



io,«-i 
t 



Q^*+l^iT,*-l 
-OX^ XjXj 



„i-l 



™«+lT,i™«-l 
X^ XjX^ 
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Rule 

I44I 



Hb\ 



Hq\ 



1471 



1481 
I49I 



Unl 



50^ 



I51I 



I52I 



I53I 



I54I 



iKcl 



55^ 



UrI 



56^ 



I57I 



I58I 



I59I 



i60i 



I6II 



l«ol 



62i 



iRQl 



63^ 



i64i 



1,02 


,«3, 


(1,1 


0,1, 


(0,0 


1,1, 


(1,0 


1,1, 


(0,1 


1,1, 


(1,1 


1,1, 


(0,0 


0,0, 


(1,0 


0,0, 


(0,1 


0,0, 


(1,1 


0,0, 


(0,0 


1,0, 


(1,0 


1,0, 


(0,1 


1,0, 


(1,1 


1,0, 


(0,0 


0,1, 


(1,0 


0,1, 


(0,1 


0,1, 


(1,1 


0,1, 


(0,0 


1,1, 


(1,0 


1,1, 


(0,1 


1,1, 


(1,1 


1,1, 


(0,0 


0,0, 



04,05,06,07) 
0,1,0,0) 

0,1,0,0 
0,1,0,0 
0,1,0,0 
0,1,0,0 
1,1,0,0 
1,1,0,0 
1,1,0,0 
1,1,0,0 
1,1,0,0 
1,1,0,0 
1,1,0,0 

1,1,0,0 
1,1,0,0 
1,1,0,0 
1,1,0,0 
1,1,0,0 
1,1,0,0 
1,1,0,0 
1,1,0,0 

1,1,0,0 
0,0,1,0 



Map 



n*+l 



r.'i-+^ r^i _L ™,* + l™*^l 






j-1 






Jb j- tXj J- Jb J- ~I~ Ju j- Ju 









t 



xlx 



i-1 



t-^t 



*+l I ^*+l^«-l ^*+l^i^«-l 



1 - x^ + xl^'x 



Xj XjX^ 



X 



t+l 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



i Xj X^ -t- X^ X( -|- XjXj 



i+l _ ™J+l™i I ^i-1 _ ™i+l^«-l 






^t-^t 



t -^t-^t 



l-x\ 

^*+l _(_ ™J 9™«+l^i ^«^«-l I ™«+l^i™«~l 

1 _ ^« + l^« _ T.«-l I ^«+l^«-l 

1 Xj X( X^ -|- X( x^ 



.i+l 



*+l™i I ™*^1 ^*+l^«— 1 



- x^ -t- x^ /:x^ Xj -t- Xj 
+2x^+lxjx^-l 

- 1 _ rJ + '^rr.i _ ry.irr.i-'^ I ^« + l^J^«-l 

- 1 Xj X^ -^t-^t "T Xj x^Xj 



24^r'+ 



X 



t+1 
^t+1 
x" 

^t+1 

^i+1 

^t+1 



J+1 ^«+l^i 



X^ Xj -|- XjXj 

j— 1 J^ „*+l„*— 1 



i—1 ™*+l^i™*— 1 



j-1 



„«+l„i^«— 1 



1 _ ^* _ t'"-^ J- T-'T-J-^'-J- _L 9^«^'-J- _ 9^«^-J-^«^« 

i t/^J. tiy-f. ~I~ tt^J. tt^J- ~I~ ZjiL/ j-JU 4- ^iJbj. Jbj-iJbj. 



*+l ^«+l^i I ^*-l ^*+l^*-l 



t+1 



Xj- XI + X 



t ^ -^t 



l-x\ + x\x\ 



Xj XjXj 



T,*+l 1 ™J 9™«+l^J 

1 _ rpi + '^rpi _ ^«-l I ^«+l™«-l I ^Irp 
i JyJ^ JyJ^ J/J- ~I~ U^J^ O^J^ ~I~ •-V-I-Uj 



«-l T,i+I^i™«-1 

•i/J- Jb J-Jb . 



t -^t-^t 



X 



t+1 



4+1 
^J+1 



- -r*+l 4- -r* — 9-T-*+l'r* -I- t-*"! — -r^+l-T-*"! — -ri-r*"^-!- 

_i_^«+l^j^«-l 
-|-Xj X(Xj 



l-xj+^x^ 



„i+l 



*+l^j^*— 1 



Xj Xj X( X^X( 
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Rule 

'651 

^66^ 



i67i 



i68i 



i69i 



1701 
I71I 
1721 
I73I 



I74I 



I75I 



lv«l 



76^ 



I77I 



^781 



I79I 



iBQi 



isU 



lool 



82i 



IfiQl 



83^ 



ifi^i 



84i 



I85I 



(ao,«i,«2: 

(1,0: 

(0,1, 

(1,1, 

(0,0: 

(1,0: 

(0,1, 

(1,1, 

(0,0: 
(1,0: 



(0,1, 
(1,1, 

(0,0, 
(1,0, 

(0,1 

(1,1, 

(0,0: 
(1,0: 

(0,1, 

(1,1, 

(0,0: 
(1,0: 



03 

0,0 

0,0 
0,0 
1,0 
1,0 
1,0 
1,0 
0,1 
0,1 

0,1 
0,1 

1,1 
1,1 

1,1 
1,1 

0,0 
0,0 
0,0 
0,0 
1,0 
1,0 



04,05,06,07) 
0,0,1,0) 

0,0,1,0 
0,0,1,0 
0,0,1,0 
0,0,1,0 
0,0,1,0 
0,0,1,0 
0,0,1,0 
0,0,1,0 

0,0,1,0 
0,0,1,0 

0,0,1,0 
0,0,1,0 

0,0,1,0 
0,0,1,0 
1,0,1,0 
1,0,1,0 
1,0,1,0 
1,0,1,0 
1,0,1,0 
1,0,1,0 



Map 



n*+l 



x\ + 2xY'x, 



t 



^' 



J+l^i-l 



+ xix 



ir^i—^ 



t-^t 









.i^i-l 



.*+l 






«+l. 






J-1 



xlx 



i-1 



t^t 



1 - 4+^ + xj+^4 



J-i 



H 



+ X 



t-^t 



«+l^«-l 



J- Jb fJu 



i-1 



t-^i 



j-1 



X 



t+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



Xt + Xt 



^*+i^«-i 
j'f J'/ 



ziX^X^ -|- Xj XjXj 



X 



+ X 



i+l^i 



«^*— 1 



x;x 



Xj Xj -t- X(Xj 



9™«+l™i™«-l 
^■^t XjX^ 



n*+l 



«+l^i 



x\ + 2xJ^-'xJ 



-OXj XjXj 



xJ^-'-xJ + X 



.*+l 



i~l 






.i+1 



^' 



.*+l 



^^^-r^ 



+ 2x^x^"^- 



i Jb-I-Jb 



i-1 



t-^i 



x! + 2Xi x! + xjx 



j„«— 1 



t-^t 



.i+1 



2xJ'*' xjx 



i-1 






n*+l 



+ X 



i+l^i 



Xi X4. 



J-1 



+ X 



i+l^i-l 



+ xlx 



i^i-1 



t-^i 



9^i+l^i™i-l 
-ZX^ XjX^ 



t+1 
i 



X. 

^t+1 
^i+1 

x" 



x: + X 



i-l 

t 



i+li-1 



Xt X 



XIX 



J-1 



t-^t 



J+1 



+ X 



i+l^i 



i+1. 






J-1 



t^t 



J+1 



™«+l™i-l 



4+1 



1 - xj + x!"^^x! 



J-1 



+ x;x 



j™j-i 



t-^t 



i+l 

/ uV-j-Jb 



i-1 



t-^i 



X 



t+1 



4+1 



4+1 



4+1 



J+1 



+ X 



i-1 



2xj+ixri 



™i™«-l I ™«+l™i„«-l 
X^X^ I -^t -^t-^t 



1 — Xl + Xf 



j+1. 



„i+l™i-l 



J+1 



X^ Xj Xj 



„i+l™i-l 






^t-^t 



H-^t 



1 — X 



i-1 
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Rule 

'861 

i87i 



i88i 



1391 



I90I 



I91I 



Inol 



92i 



lOQl 



93^ 



lo^i 



94i 



I95I 



I96I 



Invl 



97i 



^981 



lool 



99^ 



Hoo^ 



ilOli 



^102^ 



IIO3I 



H06i 



(ao,«i,«2 
(0,1 



1,1 

0,0 
1,0 
0,1 

1,1 

0,0 
1,0 
0,1 

1,1 

0,0 
1,0 

0,1 

1,1 

0,0 
1,0 
0,1 

1,1 

0,0 
1,0 



03 

1,0 

1,0 
0,1 
0,1 
0,1 
0,1 

1,1 
1,1 
1,1 

1,1 

0,0 
0,0 

0,0 
0,0 
1,0 
1,0 
1,0 
1,0 
0,1 
0,1 



04,05,06, ay j 

1,0,1,0) 

1,0,1,0 
1,0,1,0 
1,0,1,0 
1,0,1,0 
1,0,1,0 
1,0,1,0 
1,0,1,0 
1,0,1,0 

1,0,1,0 
0,1,1,0 
0,1,1,0 

0,1,1,0 
0,1,1,0 
0,1,1,0 
0,1,1,0 
0,1,1,0 
0,1,1,0 
0,1,1,0 
0,1,1,0 



(0,1,0,1,0,1,1,0) 



Map 
^t+1 - 

^t+1 - 
^t+1 - 
^t+1 - 
^t+1 - 

Table S6 



n*+i 



+ 4 



X 



i+'i-^i 



J-1 



«+l^«-l 



J-1 






+2xi+'xixi'' 



J+l^i-i 






J^i-1 



J+1. 



J-1 



J+1 



j^^+l^;*-! 



+ xlx 



i-1 



t-^t 



i+1 . 



J-1 



t-^t 



1 - x! + x!+^xj 



+ 2x|xJ-^ - 2xJ+^x^x 



J+i 



+ X 



j-1 



«+l^«-l 



2xr'x: 



1 - x! + x^+^xj 



J+i 



+ x\ 

i-l 



^«+l^«-l 



+ xlx 



j-1 



t-^i 



«+l^J^*— 1 



t-^t 



J- JU-l-^Xi 



t-^t 



^«+l^j 



/ ~r •Ji'f'Juf 



•-t 



™«+l^«-l 
J+1™J^«-1 



J+1 



+ x! 



X 



«+i^j 



Xt + Xt 



j-1 



«+l^«-l 



2xr^xJ 



x^x^ 1+ 



+x!+-^x!x 



i-l 






x'.'^'xi + X 



J+i 



«+i^«-i 



J+i 



2xJ^ xjx 



i-l 



J+1 






Q^i+I^i™*-1 



J+1 



X^ + Xt 

^j+l 



j-1 



i™.«— 1 



x;x 






^tXf 



rr^ I 9^*i-J-^« I ^«+l™«-l 



J+1 



Z,X^ X^Xj 



i-l 



Xi + X 



i+1 j-1 

t x^ 



i^*-l 



XiX 



„«+l^i^«— 1 



t^t 



J+1 



+ X 



i+l^i 



xr^ + 2x: 



i+l^i-l 
Xi 



2x!+^x!x 



i-l 



t-^t 



xj + xj-' 



J+1 



2xlx^ 
J+1 



j-1 



+ Xf x; + X 



«+i^«-i 



x;x 



i-l 



J+1. 



„i-l 



t^t 



J+l^i 



Xj + X^ 
^J+1 



«+l^«-l 



+ x;x 



i^*— 1 



t-^t 



Q™«+I^i™«-1 
OX^ XjX^ 



Xt 



•'t-^t 



Xt + 2Xj a.j 

- 4xj+^xjxJ 



J-1 



«+l^«-l 



+2xr'xrH 



+2x*x*~"'" — /i^*+i™«~.*-i 



x!^ xl + X 



(cont.) 



j-1 



9^*+l^i^*-l 
ZXj XjXj 
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Rule 

H07I 



H08l 



H09i 



Hwl 



1114 



Hl3l 



Hul 



Hl5i 



^116i 



11171 
I1I8I 



III9I 



il20i 



1124 



I122I 
I123I 



I124I 
I125I 



il26i 



(ao,«i,«2: 
(1,1 

(0,0, 
(1,0, 

(0,1, 

(1,1 

(0,0, 
(1,0, 

(0,1, 

(1,1 

(0,0, 
(1,0, 
(0,1, 

(1,1 

(0,0, 
(1,0, 

(0,1, 

(1,1 

(0,0, 
(1,0, 
(0,1, 



03 

0,1 

1,1 
1,1 

1,1 
1,1 

0,0 
0,0 

0,0 
0,0 
1,0 
1,0 
1,0 

1,0 
0,1 
0,1 

0,1 
0,1 

1,1 
1,1 
1,1 



04,05,^6,07) 
0,1,1,0) 

0,1,1,0 
0,1,1,0 

0,1,1,0 
0,1,1,0 
1,1,1,0 
1,1,1,0 

1,1,1,0 
1,1,1,0 
1,1,1,0 
1,1,1,0 
1,1,1,0 

1,1,1,0 
1,1,1,0 
1,1,1,0 

1,1,1,0 
1,1,1,0 
1,1,1,0 
1,1,1,0 
1,1,1,0 



Map 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



n«+l 



«+l^« 



xl + 2xr xj + X 



J+l^i-l 



+ xix 



«^*— 1 



t-^t 






Xj -|- X( x^ 



9^i+l^i^i-l 
ZX^ XjXj 



.«+l 



+ X 



,*+l^i 



i-l , 9 i+1 J-1 I ™i„J-l 



fJJb J- Jb -l-Jb ^ 



H-^t 



xj + xj-' 



xjx 



i-l 



t-^t 



j+1 

J- Jb-l-JU 



i-l 



.«+l 



„*+l 



+ Xj- x; + X 



t-^t 



«+i^«-i 



„i+l 



2xj XjX^ 



i-l 



n*+l 






1 - xj + x^'^ix^ 

_9^i+l^i^i-l 
Z,X^ XjX^ 



J-1 



+ X 



*+l^«-l 



+ XiX 



i-l 



t-^t 



J+^ 



+ X 



j-1 






xjx 



•i-l 



t-^t 



1 — x\ + X 



«+l^i 



j + 1 

X Uj + Jb 



i-l 



t-^t 



.«+l 



X 



+ xj 

i-l 



X'i~r J- ,^2 



XiX 



i-l 



t-^t 



+ X 



«+l^«-l 



i + 1 

y. Jbj-Ji/ 



i-l 



n*+l 



«+l^i 



+ xl — x;^ x; + X 



j-i 



t-^t 






2xjxt 



i-l 



+ 



_i_„i+l„i^i-l 
-|-Xj X^Xj 

a-xjxj-i 



.«+i 



+ xlx 



i-l 



J+l^i^*-l 



t-^t 



ZX^ XjX^ 



1 - x! + xl+^xl 



J-i 



4 

J /ya t'T -L Q^ t ,ya 

fJJbJ- iLi-l-Jb 



+ X 



*+l^«-l 



+ 2x!x' 



i-l 



t^t 



i-l 



.*+l 



t-^i 



+ X 



i-l 



«+l^«-l 



X^ X 



i+1 



i-l 



t-^t 



1 T-* J- T*+1t* J- ^«T.*~1 

± X^ ^ X^ X+ ~\~ 'JbfJb^ 



n« + l 



X 



t+1 ■ 
^t+l - 

Table S6 



X 



4 
+ x| 

i-l 



4-^i 



9^«+l™i^«-l 



X'^H" J- ,^^2 



+ X 



«+l^«-l 



/ya^~r-'- /ya^ /ya^ J- 

- X^ X^X^ 

+ x\xl~^ 



i+l„i^i-l 



.*+l 



(cont.) 



,*+i. 



+ xl — Xi Xi + X 



J-1 



^Xj XjXj 






i^*-l 



XiX 



t-^i 
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Rule 



H27I 



^129^ 



^UOl 



Hsil 



^132i 



I133I 
'134 



I135I 



H36i 



H37i 



^138i 



I139I 



^141^ 

1 

2 



142 



^43^ 



il44i 



I1 /ic:l 



145i 



^146^ 



1147^ 



H48i 



(ao,«i,«2 

(1,1 
0,0 
1,0 
0,1 

1,1 

0,0 
1,0 
0,1 

1,1 

0,0 
1,0 

0,1 

1,1 

0,0 
1,0 
0,1 

1,1 

0,0 
1,0 
0,1 

1,1 

0,0 



03 

1,1 

0,0 
0,0 
0,0 
0,0 
1,0 
1,0 
1,0 

1,0 
0,1 
0,1 

0,1 
0,1 

1,1 
1,1 
1,1 
1,1 

0,0 
0,0 
0,0 

0,0 
1,0 



04,05,06,07) 

1,1,1,0) 

0,0,0,1 

0,0,0,1 

0,0,0,1 

0,0,0,1 

0,0,0,1 

0,0,0,1 

0,0,0,1 

0,0,0,1 
0,0,0,1 
0,0,0,1 

0,0,0,1 
0,0,0,1 
0,0,0,1 
0,0,0,1 
0,0,0,1 
0,0,0,1 
1,0,0,1 
1,0,0,1 
1,0,0,1 

1,0,0,1 
1,0,0,1 



Map 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 






t-^t 






.«+! 






J-1 



«+l^«-l 



+ x'r^x 



^x\x^ 



i-\ 



•i-1 



J+l^i-1 



X^ X 



^j^«-l I 9^*+l i i-1 



^t^t 



•■"t-^t 



.«+i 



j+i. 



x; + x-T x\ + X 



'^^xix^-^ 



t-^t 



X+ X4. 



XtX 



t-^t 



+ 2x\+^x\x 



.«+i 



n«+i 



X 



i-1 



+ X 



«+l^*-l 



t-^t 



J+l. 



+ Xj- xix 



i-l 



t-^t 






i-l 



JU-I- Jb-I- 



2x\x^^+ 



H-^t 



n«+l 



+ 2x^'x\x\ 



i^i-^ 



XIX 



t-^t 



J+^ 






J-1 



«+l^«-l 



+ x.^ X 



+ 



+2x!x 



i-l 



t-^t 



?+l 

XI /y Ij rf 

-I- JUfUj 



i-l 



t^t 



J-1 



«+l^i-l 



X-t X 



J+1, 



+ X, x.x 



j-1 



t-^t 



J+1 



J+1 



J+1. 



J -t- X( -^t "T XjX^ 



j-1 



+ X/ 0/+ ~p X/ JbfJu 



4 
J+1 



i-l 



i+1. 



„i-l 



X, 






+ X 



i-l 



4-^i 
«+l^*-l 



+ ^iXt 



i-l 



2,«+l2,«-l 



^i^*-l I 9™«+l™i™«-l 



^t^t 



•"t-^t 



1 _ ^«+l I ^«+l^i^*-l 



^t-^t 



J+1 



f X/ 



x: — X 



j-i 



^«+i^«-i 

x+ x + 



i-l 



+ 2xj+ixjxr^ 



4-^i 



+ XjXj + X 



J+1 



XjX 



i-l 



t-^t 



™J _ 
2^t+l - 


- ^«+i _ 
+3x^+1 


X^ Xj + Xj 

xjxr^ 


- 2x^'x^' 


~ ^t^t 


-'+ 




™J — 

2^t+l - 


= l-xj 


— Xf x^ + 2xj x^Xj 








^t+1 - 


= A^' + 


^t ~ 2x^ Xj - 


^'^' - 


X^X-f- 


+ 3xj+^ 


^t^t 
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i-l 
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Rule 



I149I 



'1501 



^ISU 



il52i 



^153^ 



'154 
I155I 



I157I 



^158^ 



I159I 
'I6OI 



II6II 



il62i 



^163^ 



^164 






'1671 
'1681 



il69i 



(ao,«i,«2: 
(1,0, 

(0,1: 



(1,1 

(0,0: 

(1,0: 

(0,1. 

(1,1 

(0,0: 
(1,0: 

(0,1. 

(1,1 

(0,0. 
(1,0. 

(0,1. 

(1,1 

(0,0: 
(1,0: 

(0,1. 

(1,1 
(0,0: 

(1,0: 



03 

1,0 
1,0 

1,0 
0,1 
0,1 

0,1 
0,1 

1,1 
1,1 
1,1 

1,1 

,0,0 
,0,0 

0,0 
0,0 
1,0 
1,0 
1,0 
1,0 
0,1 
0,1 



04,05,06, ay j 

1,0,0,1 
1,0,0,1 



1,0,0,1 
1,0,0,1 
1,0,0,1 

1,0,0,1 
1,0,0,1 
1,0,0,1 
1,0,0,1 
1,0,0,1 

1,0,0,1 
0,1,0,1 
0,1,0,1 

0,1,0,1 
0,1,0,1 
0,1,0,1 
0,1,0,1 
0,1,0,1 
0,1,0,1 
0,1,0,1 
0,1,0,1 



Map 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



^«+l^j 

X/ X + 



„i-l 



«+l„J™i— 1 



+ 2xr'x\x[ 



n*+i 



+ Xt + x^ 



i-l 



2x:+ix: 



4 2x^ Xf- 



-2xjxri 



+ Ax\ 



*+i. 



J-i 






i+l^j-l 



X. X 






H-^t 



H-^t 



.*+l 






«+l^«-l 



Xi X 



„i^i—l 



J+1. 






J-l 



t^t 



J-1 



+ 2x\x 



i-l 



-2x\x\-^ + 4xj+i 



XtX 



i-1 



t^t 



n*+l 



n*+l 



+ X 



X 



i+1 
t 



.«+l 



-XfX 



+ X 
i-1 



i+l 
t 


Xf + Xf 


1 - 2x^+1 


xri+2xj+ixjxr 


-1 


Xf Xf + 


x|xr 


-' + 


Xf XjXj 






i 
t 


2x1-^^x1 




~|~ ZiJbf JLj-JUm 


-1 




x\ 


- 4~' + 


^t^t 


-' + 


Jb J- Jb 4- ti-i 4- 






i 
t 


2xj Xf 


+xi' 


-1 _ 


2x^+1x^-1- 







i+l. 



t-^t 



-|- OXf X^X^ 



J-1 






«+l^«-l 



^i+l™i-l 






Xi X 



J+1 



+ ~|~ -^Z '*'+ -^^ 



i-1 



+ 2x1 



i+l^i-l 



+ 



-|-XjXj x^ x^x^ 



J-1 



XIX 



J-1 



t-^t 



2 + 1 

+ X. xjx 



i-1 



t-^t 



J+1 



J+1 



J+1. 



x; + x:^ xl + X 



*+i^«-i 






i+l^«-l 



X, 



J-1 



J+1 



+ Xf XfX 



i-1 



J+1 



x; — X 



i+l 



xri+2xj 

„i-l 



«+l^«-l 



„i^«-l 



t-^t 



„i+i ^i rr.'i'—^ 



dyf ~\~ Jb f iLiJbfJbf ~\~ LiJLjf Jb f Jb f 



J+1 



+ X 



«+l^«-l 



xfx* 1 



™i 


= -j+^-r^ 


^t+1 '- 


= 1 - 4'^' 




+2xjxj"i 


Table S6 


(cont.) 


80 





+ xjx! 1 



J+1. 



+ x 

i-1 



i+l 



XiX 



J-1 



t-^t 



-t-^t 



Xf + X 



J+1 



«+l. 



„i-l 



i+l„i-l 



+2xr^xr^+ 



2xr x!x 



i-1 



t^t 



Rule 

H70i 
11711 



il72i 



^nsl 



il74i 



I175I 
'1761 



I177I 
11781 



il79i 



il80i 



II8II 
'1821 



I1 sqI 



183 



II84I 



II85I 



II86I 



il87i 



'1881 
il89i 



il90i 



(ao,«i,«2 
0,1 

1,1 



0,0 
1,0 
0,1 

1,1 

0,0 
1,0 
0,1 

1,1 

0,0 
1,0 
0,1 

1,1 

0,0 
1,0 
0,1 

1,1 

0,0 
1,0 
0,1 



03 

0,1 
0,1 

1,1 
1,1 
1,1 
1,1 

0,0 
0,0 
0,0 

0,0 
1,0 
1,0 
1,0 

1,0 
0,1 
0,1 
0,1 
0,1 

1,1 
1,1 
1,1 



04,05,^6,07) 

0,1,0,1 

0,1,0,1 



0,1,0,1 
0,1,0,1 
0,1,0,1 
0,1,0,1 
1,1,0,1 
1,1,0,1 
1,1,0,1 

1,1,0,1 
1,1,0,1 
1,1,0,1 
1,1,0,1 

1,1,0,1 
1,1,0,1 
1,1,0,1 
1,1,0,1 
1,1,0,1 
1,1,0,1 
1,1,0,1 
1,1,0,1 



Map 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



J-1 



n«+l 






t-\- Xf X^ -t- Xj 



^t-^t 



„«+l„j^«— 1 



■^t -^t -^t ' -^t ^t 



.«+l 



i-l I 9 «+l t-1 I „i™«- 



^^+l^i.>.^~l 






XiX 



i-1 



t-^t 



+ xl'^^XfX 



i-1 



t-^t 



.«+l 



+ X 



«+l^«-l 



n*+l 



J+^r^i 



^\^' 



4 



i+l„j^i— 1 






H-^t 



■■^+iTri+xjxj-i 



J-1 









xjx^ 1+ 



+2xj+^xjxj-^ 



i+l^i^i—1 



1 T*^ _L 'y*' I"!- /y^ 'y*' 

1 X( ^ X^ XjX^ 



4+1 + xl - 2xj+ixj 



XiX 



t-^i 



+ 2xJ+ix|xJ" 



X Uj-I- JUm JjJ. ~]~ tiyJ- JUM ~\~ Uj-I- KLJ-t-uVf 



H-^t 






9™i™*-l I Q^«+1^J^«- 



1 



J+1. 



+ 2xJ+ixjx^ 



J+i 



J+l^i 



XI + XiX 



«^*— 1 



J-1 



+ x 



t^t 



«+l^«-l 



+ 2x!x 



i^*— 1 



t-^t 



?+l 

X I -*- /y> *■ /y> 



j-1 



J+1 



J+l^i 



X} + X 



i-1 



^i+l^«-l 



i+l^i^j-l 



+ X^ xlx 



t-^t 



1 



X( T X^X( 



J+1 



i+l^j 



r.'i-+^ r^i ^"i-—^ 






4+1 



4+1 



1 



J+1. 



^' 



-^+1t^-i + xJxJ-i 



T-*+l 4- -r* — 9t.*+1t.* _|_ a-*~l — T-^+l-T-*"! — T.i^«-l_|_ 



4-*'* 



+2xj+^xjxj~^ 



I191I 



(1,1,1,1,1,1,0,1) 



J+1 



i+l ^i ^i—'^ 



4+1 



Xt' I /ya *■ T -^ ,'ya ^ /ya 



t-^t 
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Rule 



^1921 



I193I 



H94i 



H95i 



^196i 



I197I 
il98i 



I199I 
^200^ 



^20U 



i202i 



^203^ 
^204 



^205^ 



^206^ 

1 

1 
2 



^207.1 



^208^ 



^209^ 
I2IOI 



I2III 
I212I 



^213i 



(ao,«i,«2: 
(0,0, 

(lA 



(0,1, 

(1,1, 

(0,0, 

(lA 

(0,1, 

(1,1, 

(0,0, 

(lA 

(0,1, 

(1,1, 

(0,0, 

(lA 

(0,1, 

(1,1, 

(0,0, 

(lA 

(0,1, 

(1,1, 

(0,0, 

(lA 



03 
0,0 
0,0 

0,0 
0,0 
1,0 
1,0 
1,0 
1,0 
0,1 
0,1 

0,1 
0,1 
1,1 
1,1 

1,1 
1,1 
0,0 
0,0 
0,0 
0,0 
1,0 
1,0 



04, 05,06, ay j 



0,0,1,1 



0,0,1 



0,0,1 
0,0,1 
0,0,1 
0,0,1 



0,0,1 



0,0,1 



0,0,1 



0,0,1 

0,0,1 
0,0,1 



0,0,1 



0,0,1 



0,0,1 
0,0,1 
1,0,1 
1,0,1 



1,0,1 



1,0,1 



1,0,1 



1,0,1 



Map 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4 -^t 



n«+i 



«+i. 



: 1 - xr -xl + 2xr x: 

-|-XjXj x^ x^x^ 



^r 



^'xt'- 



x', xl + X 



j-l 



.«+l 






«+l. 



i^i-1 



XiX 



t-^t 



i+1 

+ Xj- xix 



i~l 



t-^t 



xl + 2xl^'xl 



XiX 



i-1 



t^t 



+ xl'^^xix 



i-l 



t-^t 



1 - xj+1 + xj+^xj 



J-1 



xj + xj-i 



„«+l™«-l 



+ X 



,i^i—l 



«+l^«-l 



ZiXj^X^ ~\~ ZiX 



i+l^i^i— 1 



XiX 



t-^t 



,«+l^i 



+ X 



ry'^^^ry'l' _{_ ry^'i' rf^ ^ 

Xj Xj i- x^X( 



XiX 



J-1 



t-^t 



+ xj+^xjx*-^ 



J+l^i^J-1 



1 



J+1 



xJ + 2x' 



«+l^i 



J-1 



+ X 



t-^t 



«+l„«-l 



+ 2x^x^1- 



J+1 



2x^^ xjx 



i-l 



x!^"'"x! + X 



J+1 



j-1 






J+1 



Q^'i _|_ Qq^"^"!" J- ,^.2 _1_ ,y.?' ^j-.*' 



J+1 

J. dj^Jb 



i-1 



t-^t 



J+1 



+ X 



,«+l^i 



^r^ 



„«+i_«-i 



+ xix 



i^i-l 



t-^i 



y.i + '^rr.irf.i—^ 



Xt + Xt 



j-1 _ ^j+l^j-1 



XIX 



J-1 



l^j^«-l 



t-^t 



+ X. xlx 



t-^t 



X 



,«+l^i 



X 



i+1 
t 



+ X 



«+l^«-l 



X^ X 



«+l^j^«— 1 



+ Xi xlx 



1 - xj + x!+^xj 



J-1 



J+1 



+ X 



j-1 



2xj+ixr^ 



+ x\x\'^ 



XIX 



t-^t 



i+l^i^i—1 



+ 2xr^x\xl 



T — ^i M ™«+l™i _ ^«+l^«-l I ^«+l^J™i-l 

X tt^+ ~p dy-j- Jbj^ JjJ. ^X ~y JUj^ JbfJU^ 



^«+l I ^j _ ^«+l^j _ ™«+l™«-l 

ti'J^ ~I~ O/J^ Jbf JUf dy-l- tiy-l- 



1 _ -r*~l -I- -T-^+^T-^-r^^l 
± X( "T -i-t -^t-^i 



4-^i 
J-1 



+ 2x^+^x1x^1 
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Rule 



I214I 



If 



'2 
^216i 



I217I 



^218i 



I219I 
^220^ 



I221I 
^222^ 



^223i 



^224 
^225^ 



^226^ 
^227^ 



^228i 



^229i 



^230i 



1234 
^232^ 



^233^ 



(ao,«i,«2 
(0,1 



1,1 

0,0 
1,0 
0,1 

1,1 

0,0 
1,0 
0,1 

1,1 

0,0 
1,0 

0,1 

1,1 

0,0 
1,0 
0,1 

1,1 

0,0 
1,0 



03 

1,0 

1,0 
0,1 
0,1 
0,1 
0,1 

1,1 
1,1 
1,1 

1,1 

0,0 
0,0 

0,0 
0,0 
1,0 
1,0 
1,0 
1,0 
0,1 
0,1 



04,05,06, ay j 



1,0,1,1) 



1,0,1,1 



1,0,1 
1,0,1 
1,0,1 
1,0,1 



1,0,1 



1,0,1 



1,0,1 



1,0,1 
0,1,1 
0,1,1 



0,1,1 



0,1,1 



0,1,1 
0,1,1 
0,1,1 
0,1,1 



0,1,1 



0,1,1 



Map 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



4+1 



n*+l 



+ 4 



xl^^xl + X 



i-1 



n«+l^«-l 



i-1 



2xl^'xr' - 2x1x1'' + 



-1-OX^ Xj-X^ 



J+l^i-i 






J^i-i 



J+1 



ZX^ XfX^ 



.«+l 



x^^^x^'-^ 



l-xi + xi+^xi 



+ XtXt 
i-1 



j-1 



X 



+ 2x1x1 



i-l 






t-^t 



.*+l 



+ X 



i-l 



J+l^«-l 






«+l. 



„i-l 



t-^t 



I- xi + xl~^^xi 



n*+l 



4 
+ xl 

i-l 






+ xlx 



i-l 



t-^t 



Jb-t Jb + 






j+l„j„« — 1 



+ x-f xix 



t-^t 



Xf -|- X^Xf- 



n*+l 



+ xl 



X 



i+'i-^i 



i-l 






X^ X 



Xi + Xt 



i+l^j^i— 1 



«+l^«-l 



2x^^x1 



44 ^+ 



+ Xl^ XfX 



x'.'^^xi + X 



.«+i 



«+i^«-i 



t-^t 



«+i^j^*— 1 



-I- tiy-j-Uj 



t-^t 



.i+1 






-z,x^ X^X^ 



n* + l 



x^ + x^ 



j-1 



«„.«—! 



x;x 



J+'i-^i 






J+1. 



J-1 



Xi + X 



j+1 i-l 
t ^t 



«^*— 1 



XiX 



t^t 



.«+l 



+ X 



«+l^i 



J-1 



+ 2xr x: 



i+l^«-l 






t-^t 



xl + xl' 






4-^t 



4-^4 



J+1 



+ x!"^-^xj + X 



«+i^«-i 



x;x 



i-l 



t-^t 



J+l^i 



1 



Xj + X^ 
^i+1 



«+l^«-l 



+ x;x 



j^*— 1 



t-^t 



Z,X^ X^X^ 
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Figure S16: Spatiotemporal evolution of the 256 Wolfram rules R2 [R indicated over each 
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Figure S16 (cont.): Spatiotemporal evolution of the 256 Wolfram rules R^ 
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Figure S16 (cont.): Spatiotemporal evolution of the 256 Wolfram rules R^ 
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Figure S16 (cont.): Spatiotemporal evolution of the 256 Wolfram rules R^ 
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Figure S16 (cont.): Spatiotemporal evolution of the 256 Wolfram rules R^ 
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Figure S16 (cont.): Spatiotemporal evolution of the 256 Wolfram rules R^ 
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Figure S16 (cont.): Spatiotemporal evolution of the 256 Wolfram rules R^ 
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Figure S16 (cont.): Spatiotemporal evolution of the 256 Wolfram rules i?2 
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